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Abstract 

A comprehensive review of the discrete quantum mechanics with the pure imagi- 
nary shifts and the real shifts is presented in paraUel with the corresponding resuhs in 
the ordinary quantum mechanics. The main subjects to be covered are the factorised 
Hamihonians, the general structure of the solution spaces of the Schrodinger equation 
(Crum's theorem and its modification), the shape invariance, the exact solvability in the 
Schrodinger picture as well as in the Heisenberg picture, the creation/annihilation oper- 
ators and the dynamical symmetry algebras, the unified theory of exact and quasi-exact 
solvability based on the sinusoidal coordinates, the infinite families of new orthogonal 
(the exceptional) polynomials. Two new infinite families of orthogonal polynomials, 
the Xi Meixner-Pollaczek and the Xi Meixner polynomials are reported. 

PACS: 03.65.-w, 03.65. Ca, 03.65.Fd, 03.65.Ge, 02.30.Ik, 02.30.Gp 

keywords: discrete quantum mechanics, difference Schrodinger equation, exact solv- 
ability, quasi-exact solvability, intertwining relations, shape invariance, Heisenberg operator 
solutions, closure relations, new orthogonal polynomials 



1 Introduction 

In the discrete quantum mechanics (dQM), developed by the present authors, the time- 
independent Schrodinger equation is a second order difference equation instead of a second 
order differential equation in the ordinary quantum mechanics (oQM). The basic framework 



of quantum mechanics, such as the probabihty interpretation, the Hilbert space structure, 
the continuous real time, etc. is unchanged and the only difference is the discrete forms of 
various Hamiltonians. The discrete quantum mechanics can be considered as a generalisation 
(deformation) of the ordinary QM in the sense that a differential equation can be obtained 
from a difference equation in an appropriate limit. The generalisation is also consistent with 
Heisenberg's idea that a certain fundamental length would appear in the theory at some stage 
of descending the microscopic ladder. This generalisation has turned out to be extremely 
fruitful. Most of the concepts and methods of oQM can be easily transplanted to dQM, of- 
fering a unified platform for the systematic understanding of various orthogonal polynomials 
satisfying second order difference equations with the pure imaginary shifts (idQM) as well as 
the real shifts (rdQM). Sometimes these polynomials are said to have the bi-spectral property 
[1]. They satisfy two types of spectral conditions, the difference eigenvalue (Schrodinger) 
equation (12.231) and the three term recurrence relations (I2.139p . Those polynomials satisfy- 
ing the difference equations with the real shifts are also called the orthogonal polynomials of 
a discrete variable [2]. They contain most of the (g)-hypergeometric orthogonal polynomials 
belonging to the Askey scheme [21 [3l 01 El |6] . The power of the discrete QM is demonstrated 
most eloquently in the discovery of various infinite families of new orthogonal (exceptional) 
polynomials. In less than two years after the discovery of the infinite family of the exceptional 
Laguerre and Jacobi polynomials [71[8l[9], the generic counterparts in dQM, the exceptional 
Wilson, Askey- Wilson, Racah and g-Racah polynomials are constructed by the present au- 
thors flU\ [TT] . These new orthogonal polynomials are expected to play a significant role in 
mathematics, mathematical physics and related disciplines. The exceptional Jacobi polyno- 
mials provide infinitely many global solutions of Fuchsian differential equations with more 
than three and arbitrarily many regular singularities |2] , although the locations of the extra 
singularities are rigidly specified. 

This topical review of dQM provides a comprehensive overview of the subjects related to 
the exactly and quasi-exactly solvable quantum particle dynamics [I2]^|2S]- Due to the length 
constraints, we have to concentrate on the systems of single degree of freedom, which is the 
most basic and best established part of the theory. The subjects to be covered in this article 
are the adaptation to dQM of the methods and concepts of oQM accumulated over 80 years 
after the birth of quantum mechanics. After the general setting of the discrete quantum 
mechanics with the pure imaginary shifts (idQM) as well as the real shifts (rdQM), we start 
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with the factorised Hamiltonians and the Schrodinger equations [121 ISQl |2T]. The general 
structure of the solution spaces is explored by the intertwining relations and Crum's theo- 
rem together with its modification [27j-[29], [21]^[26]. Exact solvability in the Schrodinger 
picture is explained by the shape invariance [301 Ell [121 |20l [21]. The generic eigenvalue 
formula, unified Rodrigues formulas and the forward/backward shift operators are deduced. 
The solvability in the Heisenberg picture is derived based on the closure relation between the 
sinusoidal coordinate and the Hamiltonian [131 [T5] . The creation/annihilation operators are 
introduced and their connection with the three term recurrence relations of the orthogonal 
polynomials is emphasised. The dynamical symmetry algebra generated by the Hamiltonian 
and the creation/annihilation operators are also established for all the solvable systems in 
the Heisenberg picture [2n [23]. This includes the dynamical realisation of the g-oscillator 
algebras [19]. For the rdQM, the dual pair of orthogonal polynomials (the Leonard pair 
[321 [33] ) and the related Askey- Wilson algebras [,34j are explored. The unified theory of 
exactly and quasi-exactly solvable dQM is another major theme of this review. Based on 
the fundamental properties of various sinusoidal coordinates, a simple recipe to construct a 
solvable dQM Hamiltonian is provided [23]. It generates all known exactly solvable theories 
and many new ones. As a byproduct it also gives a recipe of constructing quasi-exactly 
solvable dQM Hamiltonians. The final major subject of the review is the new (exceptional) 
orthogonal polynomials [7] [TT] . [35] [38] . 

Throughout the review we stress the similarity and differences among oQM, idQM and 
rdQM. As far as possible we adopt the same notation and present one all-inclusive formula 
for the three categories, oQM, idQM and rdQM. The differences are exhibited by explicit 
examples. We usually adopt three examples from each group denoted by the name of the 
corresponding eigenpolynomials with an abbreviation. For the oQM, they are the Hermite 
(H), the Laguerre (L) and the Jacobi (J) polynomials. For the idQM, we choose the special 
case of the Meixner-PoUaczek (MP) (the parameter (p is fixed to (p = |), the Wilson (W) and 
the Askey- Wilson (AW) polynomials. For the rdQM, the Meixner (M), the Racah (R) and 
the g- Racah (gR) polynomials are selected. In some cases, we present only one representative 
from each group. They are usually the most generic ones, J, AW and gR or the simplest ones 
H, MP and M. For the full line-ups, the original references should be consulted [201 [211 [23]- 
We will provide cross references to the original articles with the section or the equation 
numbers. 
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This paper is organised as follows. In section two various concepts and methods in dQM 
are explained. These are: (a) the factorised Hamiltonians (I2.5p - fl2.12p & the 'Schrodinger 
equations' for the eigenpolynomials fl2.23p - fl2.27p together with the groundstate eigenfunc- 
tions f l2.52p . f l2.55p as the orthogonality weight functions fl2.49p -f l2.50p . (b) the intertwining 
relations connecting various Hamiltonians fl2.58p . fl2.65p -f l2.66p . and eigenfunctions f l2.6ip . 
f l2.68p . (12. 73 p . This is a brief summary of Crum's theory in our language. Its deformation 
a la Krein-Adler is also mentioned, (c) shape invariance f l2.105p and the full energy spec- 
trum fl2.106p and the unified Rodrigues type formula fl2.107p for all the eigenfunctions. (d) 
solvability in the Heisenberg picture fl2.144p and the closure relation fl2.142p . The formulas 
for the creation/annihilation operators fl2.147p - fl2.148p are given, (e) dual closure relation 
f l2.182p together with its connection with the Askey- Wilson algebra. In section three we 
provide the essence of the unified theory of exact and quasi-exact solvability in dQM. In sec- 
tion four the new (exceptional) orthogonal polynomials satisfying second order differential 
(difference) equations are explored. The exceptional Meixner-PoUaczek and the exceptional 
Meixner polynomials are new results. The final section is for a summary and comments. 
Basic symbols and definitions are listed in Appendix. A substantially shorter version of the 
present review is published recently [39] • We apologise to all the authors whose good works 
could not be referred to in the review due to the lack of space. 

2 Discrete Quantum Mechanics 

The dynamical variables of one-dimensional QM are the coordinate x and its conjugate 
momentum p, which is realised as a differential operator p = —i^^ = —ihdx- Hereafter we 
adopt the convention h= 1. The dQM is a generalisation of oQM in which the Schrodinger 
equation is a difference equation instead of differential in the oQM [I2]^[I1], [201 EI]- In 
other words, the Hamiltonian contains the momentum operator in exponentiated forms e^^^ 
which work as shift operators on the wavefunction 

For the two choices of the parameter /3, either real or pure imaginary, we have two types of 
dQM; with (i) pure imaginary shifts = 7 G M^o (idQM), or (ii) real shifts (5 = i (rdQM), 
respectively. In the case of idQM, '0(x =F Z7), we require the wavefunction and potential 
functions etc to be analytic functions of x with their domains including the real axis or 
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a part of it on which the dynamical variable x is defined. In contrast, in the rdQM, the 
difference equation gives constraints on wavefunctions only on equally spaced lattice points. 
Then we choose, after proper rescaling, the variable x to take value of non- negative integers, 
with the total number either finite (A^ + 1) or infinite (xmax = N or oo). To sum up, the 
dynamical variable x of the one dimensional dQM takes continuous or discrete values: 

idQM: xeR, xe{xi,X2); rdQM: x G Z>o, 2;G[0,Xmax]. (2.2) 

Here Xi, x^ may be finite, —oo or +00. Correspondingly, the inner product of the wavefunc- 
tions has the following form: 

idQM: U.9)= / f{xyg{x)dx- rdQM: (/, f?) = /(a;)*(/(x), (2.3) 

and the norm of f{x) is ||/|| = a/(/, /). 

We will consider the Hamiltonians having a finite (rdQM with finite A^) or semi-infinite 
number of discrete energy levels only: 

= ^(0) <^(1) <£(2) < ■■■ . (2.4) 

The additive constant of the Hamiltonian is so chosen that the groundstate energy vanishes. 
That is, the Hamiltonian is positive semi- definite. It is a well known theorem in linear algebra 
that any positive semi-definite hermitian matrix can be factorised as a product of a certain 
matrix, say A, and its hermitian conjugate . As we will see shortly, the Hamiltonians we 
consider always have factorised forms in one-dimension as well as in higher dimensions. 

2.1 Factorised Hamiltonian 

The Hamiltonian we consider has a simple factorised form [10] 

D 

'H = A^A or n = ^A]Aj in D dimensions. (2.5) 



The operators A and A^ in one dimensional QM are: 

oQM: ^,^d_dMx^^ ^ 0o(x) = e-(^), (2.6) 

dx dx dx dx 

'H=p^ + U{x), Uix) = {d,w{x)y + dlw{x), (2.7) 
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idQM: A = ^{e^P^/V*{x) - e-iPy/V{x)), 7 e M^o, 

= -i{^/V{x)eiP ~ ^/W{x)e-^P), V{x),V*{x) e C, (2.8) 

n = ^/V{x)e^P^/V*{x) + y/V*{x)e-^P^/V{x) - V{x) - V*{x), (2.9) 

rdQM: A= y^B{x)-e^^/D{x), A^ = ^/B{x) - ^/D{x) , (2.10) 

D{x) > (for a; > 0), D(0) = 0, 

S(x) > (for X > 0), 5(A^) = (for the finite case), (2.11) 

n = -^B{x)e^^D{x) - ^yD{x)e-^^/B{x) + B{x) + D{x). (2.12) 

Here w{x) is called prepotential. The function V*{x) in idQM is an analytic function of 
x obtained from V{x) by the ^-operation, which is defined as follows. If f{x) = ^ 

n 

an G C, then f*{x) =^ ^a*x", in which a* is the complex conjugation of a„. Obviously 

n 

f**{x) = f{x) and f{x)* = f*{x*). If a function satisfies /* = /, then it takes real values on 
the real line. To the best of our knowledge, the factorised Hamiltonian approach to dQM 
is our invention. The existing factorisation methods for dQM are not for the Hamiltonian 
"H, but for Hamiltonian minus an eigenvalue "H — X{n) H2] . The condition D{0) = in 
rdQM (12.111) is necessary for the term '0(— 1) not to appear in 7^-0(0). Likewise B{N) = is 
necessary for the finite case. Roughly speaking these correspond to the regular singularities 
in oQM. The Hamiltonians of idQM and rdQM can be written in a unified notation: 



dQM: n = e(^^/V+{x) e^P^/V4x) +^/V4x) e-^^^/V+ix) - V+{x) - V^-(x)j, (2.13) 
idQM: /3 = 7, £ = 1, V+{x) = V{x), V^{x) = V*{x), 
rdQM: /3 = z, £ = -1, V+{x) = B{x), V^{x) = D{x). (2.14) 

Multi-particle exactly solvable systems can be constructed in a similar way. The prepo- 
tential approach is also useful in Calogero-Moser systems |13] in oQM H5]: 



We also adopt the factorised Hamiltonian for the discrete analogues of the Calogero-Moser 
systems: 

idQM : Aj = z(e^^'^ ^F;(x) - e-^^'^yV^) , V,{x), V*{x) G C, 

= _ (^• = l,...,D), (2.16) 



which is markedly different from the approach of Ruijsenaars et al [16]. The multiparticle 
version of rdQM is now being developed. 
The Schrodinger equation 

H(t)n{x)=£{n)(t>n{x) {n = Q,l,2,...), = £(0) <^(1) <£(2) < ■■■ , (2.17) 

is a second order differential (oQM) or difference equation (dQM) and the groundstate wave- 
function 00 (2^) is determined as a zero mode of the operator A {Aj) which is a first order 
equation: 

A(j)o{x) = {AM^) = 0, J = 1, . . . , D) ^ HUx) = 0. (2.18) 
By definition all the eigenfunctions are required to have a finite positive norm: 

< (0„,0„) < cx). (2.19) 

Throughout this paper we look for the eigenfunctions in a factorised form: 

cj)n{x) = (t>o{x)Pn{'n{x)), (2.20) 

in which Pn{r]{x)) is a degree n (except for those discussed in ^ polynomial in the sinusoidal 
coordinate r]{x). However, some basic results, such as the (modified) Crum's theorem in 
§ I2.2H2.3] and the shape invariance in § l2.4l hold without this assumption f l2.20p . The explicit 
forms of the eigenfunctions will be given for various examples f l2.30p - fl2.39p . Here we will 
briefly introduce their general property for dQM. The sinusoidal coordinate 77(2;) in dQM is 
a monotone increasing function of x with the initial (boundary) condition 

dQM : 7^(0) = 0. (2.21) 

With this, we also normalise the eigenpolynomials {Pn{r])} in rdQM as: 

rdQM: F„(0) = 1 (n = 0,l,...). (2.22) 

In the multi-particle case n will be a multi-index. Here we require (l)o{x) to be chosen real 
and positive for the physical values of x, which is always guaranteed in oQM. One important 
distinction between a differential and a difference equation is the uniqueness of the solution. 
For a linear differential equation, the solution is unique when the initial conditions are 
specified. In contrast, a solution of a difference Schrodinger equation f l2.17p multiplied by 
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any periodic function of a period i(3 (27 for idQM, 1 for rdQM) is another solution. In the 
case of idQM, this non-uniqueness is removed by requiring the finite norm condition for the 
groundstate wavefunction (0O)0o) < 00 and the hermiticity (self-adjointness) condition 
of the Hamihonian. For details, see |21j Appendix A. For the rdQM, the periodic ambiguity 
of period 1 is harmless since the values of eigenf unctions on the integer lattice points only 
count for the inner products, etc. It is also important to realise that the periodic ambiguity 
does not appear in the polynomial solutions of fl2.23p . 

It is important to stress that the second order equation for Pn{r]{x)) is square root free: 

nPMx)) = £in)PM^)). (2.23) 

In other words, the similarity transformed Hamiltonian T-L in terms of the groundstate wave- 
function (j)o{x) has a much simpler form than the original Hamiltonian "H: 

?^ = 0o(a;)-'o-Ho0o(x), (2.24) 
oQM: H = -— -2— ^ — , (2.25) 

dQM: n = e(v+{x){e'^P-l) + V4x){e-^P-l)^ (2.26) 

V{x){e"'P-l) + V*{x){e-"'P-l) : idQM 
B{x){l - e^) + D{x){l - e-^) : rdQM ' 



(2.27) 



For all the examples discussed in this section, T-L is lower triangular 

'Hr]{x)'^ = £{n)rj{xY + lower orders in rj{x), (2.28) 

in the special basis 

1, ri{x), r]{xf, ...,77(x)", 
spanned by the sinusoidal coordinate r]{x). This situation is expressed as 

■HVn C v., V„ = Span [l, 77(2;), . . . , t^{xY] . (2.29) 

Here are some explicit examples. For the oQM the prepotential w{x) determines the 
potential U{x) of the Hamiltonian "H = + U{x), U{x) =^ [dxw{x)Y + d1w{x): 

H : w{x) = —\x^, —00 < a; < 00, 

U{x) = - 1, 7]{x) = X, (2.30) 
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L : w{x) = —\x'^ + g\ogx, 9 > 0, < x < oo, 

U{x)^x^ + ^^^^-l-2g, rj{x)^x\ (2.31) 

X'^ 

J: w{x) — glogsmx -\-h\ogcosx, g > 0, h > 0, < a; < |, 

U(x) = + _ . ^ .X ^ 2x. (2.32) 

sin x cos^x 

Let us note that a; = for L and x = 0, | for J are the regular singular points of the Puchsian 
differential equations. The monodromy at the regular singular point is determined by the 
characteristic exponent p: 

Mp = e^""'". (2.33) 

The corresponding exponents are p — g, 1 — 51 for L and p — g,l — g and p — h,l — h ior J. 
For the idQM {0 < q < 1): 

MP: V{x)^a + ix, a > 0, -00 < x < 00, 7=1, ri{x) ^ x, (2.34) 

n^=i(% + ix) 

2ix{2ix + 1) ' 

rj{x)—x'^, {al,al,al,al} — {ai,a2,a3,a4} (as a set), (2.35) 
(1 - a e'^) 

AW: V(x)^-- — , ,^ laJ < 1, < x < tt, 7 = logg, 

77(2;) = 1 — cos X, {a^, 03, ttg, 04} = {oi, 02, 03, 04} (as a set). (2.36) 



W : V{x) = I,^, , , , Re{aj) > 0, < x < 00, 7 = 1, 



For the rdQM (0 < g < 1): 

c 1 

M: B{x) = (x + /5), D(x) = x, /3 > 0, < c < 1, 

1 — c 1 — c 

r/(a;) = x, Xmax = 00, (2.37) 
X + a){x + b) {x + c) (.T + d) 



R : S(a;) = 



{2x + d){2x+ 1 + d) 



{x + d- a){x + d- b){x + d - c)x jdef , , , , , 
= (2.-l + c/)(2. + d) ' d^- + b + c-d-l, 

a ^ -N, a + b > d > 0, < c < 1 + d, 

r]{x) = x{x + d), Xmax = N, (2.38) 
(1 — aq^){l — bq^){l — cq^){l — dq^) 



qR : B{x)= „ ^ ^ .^^^^^ 



{1 - dq^''){l - dq^ 

a-Mg^)(l - fe-Mg^)(l - c-^dq^){l - q^) jdef 
(l-dg2x-i)(l-(ig2^) 
a = g"''^, < a6 < d < 1, gd < c < 1, 
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r]{x) = {q-^ -l){l-dq^), x^^ = N. (2.39) 

The R and qR. cases admit four different parametrisations specified by e, e' = ±1, see pO] . 
Here we present e = e' = 1 case only. It should be emphasised that the sinusoidal coordinates 
T]{x) in rdQM in general depends on the parameters in contradistinction to those of the oQM 
or idQM. 

The eigenf unctions of these examples have the form (12.201) . The eigenpolynomials Pn{ri{x)) 
are the following (for notation, see [6]): 
oQM : 



H : Pnivix)) = H„,ix) = (2x)"2Fo 



n n— 1 



L: PM-)) = Lt'-\x^)'='^-^±p^^FJ-';, 

nl ^9+2 



x" 



X 



J : P„(r,(x)) = ri»--'-4'(cos 2,r) S <»±i^,F, r"' " + " + " 



9 + \ 



1 — cos 2x 



idQM 



MP: PM^)) = PlC\x■l)''^''-^^\F, 



{2a)n (-n, a + ix 



nl 



2a 



W: Pn{r]{x)) = Wn{x^ ■,ai,a2,as,a4) (bi =^ ai + 02 + as + 04) 



def 



fli + a2)n(ai + a3)n(ai + a4)n 



dcf 



AW: P„{ri{x)) = Pn{cosx;ai,a2,a^,a4\q) (64 = 01020304) 

'q~"', b4q"-~^, aie^^, Oie 
O1O2, OXO3, 0^04 
rdQM : 



-n, n + 61 — 1, Oi + ix, Oi — ix 

Oi + O2, Oi + O3, Oi + O4 
def 



Oi "(O1O2, O1O3, O1O4 ; q)n 403 



M: F„(77(x)) = M„(x;/3,c)'=^%Fi( ""j^ | 1 - 
R: Pn{i]{x)) = Rn{x{x + d) ; a — 1, d — a,c — 1, d — c) ^ 4F3( ' ~'~ 'r ' 



(2.40) 
(2.41) 

, (2.42) 

(2.43) 
l),(2.44) 

(2.45) 
(2.46) 



1 

(2.47) 



qR: Pn{r]{x)) = Rn{q '' + dq^;aq ^,da ^,cq ^,dc ^\q) = ihC^ ^ ^q , q , dq 

V a, 0, c J 

(2.48) 

Obviously, the square of the groundstate wavefunction 4>o{xY provides the positive defi- 
nite orthogonality weight function for the polynomials: 



oQM, idQM: / Mx?Pn{vix))Pn{vix))dx = hJn^, 



(2.49) 



Xl 
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^max -I 

rdQM : Yl U^?Pn{'n{x))PMx)) = - 



x=0 



Sr, 



(2.50) 



Let us emphasise that the weight function, or (l)o{x) is determined as a solution of a first 
order differential (difference) equation fl2.18p . without recourse to a moment problem. This 
situation becomes crucially important when various deformations of orthogonal polynomials 
are considered. For the oQM, the weight function is simply given by the prepotential (poix)"^ = 
q'^w{x) ^YiQ normalisation constants are 

V{n + g+\)V{n + h + \) 



:H 
: L 



(2.51) 



2n\{2n + g + h)T{n + g + h) 

The explicit forms of the squared groundstate wavefunction (weight function) 0o(3^)^ and the 
normalisation constants /i„ for the above examples in pure imaginary shifts dQM are: 



Via + ix)T(a — ix) 



00 (a;) 



{ttj + ix)r{aj — ix) 



: MP 
: W 



(2.52) 



(e=^*^g)oo(e-''^g)oon(Ke'^?)oo(«,e-^^g)oo) ' : AW 
i=i 

( 27r(2=^"n!)"^r(r2 + 2a) : MP 

2?™! {n -\- bi — 1 

)n n + tti + aj) ■ r{2n + bi)-^ :W 

l<i<j<4 

27r(M""';?)n(M'";g)oo(g"+';g);o' H («^%?";?)- ^ AW 

l<t<j<4 



(2.53) 



For the rdQM, the zero- mode equation A<f>o{x) = f l2.18p is a two term recurrence relation, 
which can be solved elementarily by using the boundary condition 02.111) : 



y=o 



The explicit forms of 0o(^)^ ^"^^ d"^! which are related by duality, are: 

r (/3).c- 

{a,b,c,d)x 2x + d 



: M 



xi 



[1 + d — a,l + d — b,l + d — c, 1)^ d 
{a,b,c,d;q)x 1 - dq'^'' 

[a^^dq,b~^dq,c~^dq,q;q)xd^ l — d 
11 



: R 



: qR 



(2.55) 



(/3)nC" 



nl 



X (1 - cf 
(a, b, c, d)j 



2n + d 



[1 + d — a,l + d — b,l + d — c, 1)„ d 
(-1)^(1 + rf-a,l + rf-f',l + c^-c; 



{d+l)N{d+l)2N 

(a, b,c,d; q)n 1 — c?(?^"' 



M 



R 



(2.56) 



[a~^dq,b~^dq,c~^dq,q;q)nd^ 1 — d 



-1 



-Mg, b-^dq, c-'dq ; g);v rf^g^^^^+i^ 



gR 



{dq;q)N{dq;q)2N 

It should be stressed that ^0(2;)^ for the rdQM (12.541) . when continued analytically in x, 
vanish on the integer points outside the defining domain, i.e., x G Z<o and x G '^>xn^^^- 

To the best of our knowledge, the first example of a factorised idQM Hamiltonian was 
mentioned in 2001 eq.(4.28) of [47j for the Meixner-Pollaczek polynomial. But it was not 
considered general principle. 

2.2 Intertwining Relations: Crum's Theorem 

The general structure of the intertwining relations works equally well for the oQM as well as 
the dQM. Let us denote by "Ht"' the original factorised Hamiltonian and by "H'^^ its partner 
(associated) Hamiltonian obtained by changing the order of A'^ and A: 



dcf 



A^A, AA^ 



(2.57) 



One simple and most important consequence of the factorised Hamiltonian (12. 5p is the in- 
tertwining relations: 



= AA^A = -Hl^U, AW^ = A^AA^ = H[°U^ 



(2.58) 



which are equally valid in the oQM and the dQM. The pair of Hamiltonians 7^'°^ and "H^^' 
are essentially iso-spectral and their eigenfunctions {(j)n\x)} and {0n^(a;)} are related by the 
Darboux-Crum transformations j48[ [27] : 

(2.59) 
(2.60) 

(2.61) 

(2.62) 



n^'^^^\x)=£{n)<P^^\x) (72 = 0,1,...), ^0S"(x)=O, 
n^'^^J:\x)=£in)<P^^\x) (72=1,2,...), 



A. 

Sin] 



0W(x) (n=l,2,...), 



(0W, 0W) = S{n)i<Prr, <Pm) [u, m = 1, 2, . . .). 
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The partner Hamiltonian "H'^l has the lowest eigenvalue £{!)■ If the groundstate energy £(1) 
is subtracted from the partner Hamiltonian H'^^, it is again positive semi-definite and can 
be factorised in terms of new operators A^^'^ and ^t^'^: 

T^W = ^Wt^[i] + £(1)^ = 0. (2.63) 

It should be stressed that in rdQM with a finite A^, the size of the Hamiltonian decreases by 
one, since the lowest eigenstate is removed, see [26] . 

By changing the orders of A^^^"^ and A^^\ a new Hamiltonian 7{[^1 is defined: 

7/[2]t^f^[iUWt + £:(!). (2.64) 

These two Hamiltonians are intertwined by A^^^ and 

_ ^(1)) = ^W^WUW = (-Hp] _ ^(l))^W, (2.65) 
^[i]t(^[2] _ ^(1)) = ^[ilUHUWt = _ £(i))^Wt. (2.66) 

The iso-spectrality of the two Hamiltonians H^^^ and T/t^' and the relationship among their 
eigenfunctions follow as before: 

H[2l0L'k^)=^H0k''(^) (n = 2,3,...), (2.67) 

<P^^K^)=A^'¥^\x), 0M(x)= f"^U, 0L^k^) (r^ = 2,3,...), (2.68) 

£(n) -£(1) 

(0l?UI^^) = (^(^) - ^(1)) (0L'UW) (n, m = 2, 3, . . .), (2.69) 

^[2] ^ ^[2]t^[2] ^ ^(2)^ ^[2]0p](a;) = 0. (2.70) 

This process can go on indefinitely by successively deleting the lowest lying energy level: 

def ^[.-i]^[s-i]t + ^(3 _ 1) ^ ^WUH + S{s), (2.71) 

?/H0l:^(^)=^^(^)0l:^(^) (n = 5,. + l,...), ^W0i^'(^)=O, (2.72) 

</)W(a;)=^[^-iVr'k^), 0t^'(^) = , 7T0[:'(^) (n = s,s + l,...), (2.73) 

t[n) — t[s — 1) 

(0H, ^l:') = {£{n) - £{s - 1)) {ct>t'\^t'^) (n, m = . + 1, . . .). (2.74) 
The quantities in the s-th step are defined by those in the (s — l)-st step: (s > 1) 

oQM: u;W(x) = log|0W(x)|, (2.75) 
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idQM : V^' 

A^' 

rdQM : B^' 

A^' 



del 



x] 



dcf 



+ l)D[*-i](a; + 1) 



x + l) 



(2.76) 
(2.77) 

(2.78) 
(2.79) 

(2.80) 



^SW(x) -e^-yDH(x), = ,y^^H(^e-9. (2.8I) 



The eigenfunctions at the s-th step have succinct determinant forms in terms of the Wron 
skian (oQM) and the Casoratian (dQM) [JTlllHlEllEg: (n>s>0) 



oQM: W[A,...,/^](x)^=^det(^^^) 



\x 



dx^ / l<j,k<m 

W [00, (;/>!, . . . ,(l)s-l,(l)n]ix) 

W [00,01, ••• ,</>s-i](a^) 
idQM : W^[/i, . . . , UKx) = z^™(™-^) det (/fc(x + z 



(Wronskian), 



l<j,k<m 



(2.82) 
(2.83) 

(Casoratian), (2.84) 



j=Q 



W^[0O,01,...,0g_l,0„](x) 
W^[0O,01, • • . ,0s-l](x - i|)' 



s-1 



[X 



rdQM : Wc[/i, . . . , /J(x) = det ('/fc(x + j - 1)) (Casoratian), 

Wc[0o,0i,...,0^-l,0n](x) 
Wc[0O,01,...,0s-l](a;+ 1) 

Wc[0o,0i,...,0^_l,0„](x) 

Wc[0o,0i, . . . ,0s-i](a;) 



fc=0 
s-1 



1)^ JJ^DW(x + s- A;) 



(2.85) 
(2.86) 
(2.87) 

(2.88) 



fc=0 



The norm of the s-th step eigenfunctions have a simple uniform expression: 



s-l 



j=0 



(2.89) 



This situation of the Crum's theorem is illustrated in Fig. 1. 



A quantum mechanical system with a factorised Hamiltonian H = A'^A together with 
the associated one "H'^^ = AA'' is sometimes called a 'supersymmetric' QM [501 EI]- This 
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appears rather a misnomer, since as we have shown the factorised form is generic and it 
imphes no extra symmetry. The iso-spectrahty is shared by all the associated Hamiltonians, 
not merely by the first two. In this connection, the transformation of mapping the s-th to the 
[s + l)-st associated Hamiltonian is sometimes called susy transformation. It is also known 
as the Darboux or Darboux-Crum transformation. Those covering multi-steps are sometimes 
called 'higher derivative' or 'nonlinear' or 'A/'-fold' susy transformations [52| l53l |5i 



^(3) 



sm 



6 

00 

n 



A 



A^ 
A 



¥-< 1 

A^ 


0W ^wt 




A 








A^ 






A 







Till] 



Figure 1: Schematic picture of the Crum's theorem 
2.3 Modified Crum's Theorem 

Crum's theorem provides a new iso-spectral Hamiltonian system by deleting successively the 
lowest lying levels from the original Hamiltonian systems H and The modification of 

Crum's theorem by Krein-Adler [281 [29] is achieved by deleting a finite number of eigenstates 
indexed by a set of non-negative distinct integer^l V =^ {di,d2, ■ ■ ■ ,de} C Z>q, satisfying 
certain conditions to be specified later. After the deletion, the new groundstate has the label 



^ Although this notation dj conflicts with the notation of the normalisation constant (i„ in (|2.50p . we 
think this does not cause any confusion because the latter appears as ^ Snm- 



15 



given by 



= mm{n \ n e Z>o\V}. (2.90) 



Corresponding to fl2.75l) - fl2.88p . the new iso-spectral Hamiltonian system is [25| 126]: 

'H = A^A + £ifi), A4>^{x) = 0, (2.91) 

nM^) = S{n)Mx) {n e Z>o\V), (2.92) 

oQM: n = p^ + U{x), U{x) = U{x)-2dl(\ogW[(Pd,Ad,,...AdM^)), (2-93) 

A'^= - dr,w{x), A^ = -da: - dj:w{x), w(x) log |0^(x)|, (2.94) 

7 / N def W[(pdi,<Pd2,---,4>di,An]ix) 

(Pn{X) — r 7 , -1/ X ) [Z.JO) 

W [0rf,,0d2,...,0rfJ(a;) 



idQM: V{x) = yjv{x - 2|7)V*(x - 2^7) 



^ [(/)rf,,...,0rfj(3: + z|) [0rf,,...,0rf,,0^](x-27) 
. . . ,0dj(x - «|) [0^1, . . . ,(/)rf^,(/)^](x) 



= -2(^^y(x)e2P- ■^y*(x) e-2P), (2.97) 



, . . . , (a; - i|) [c;/)^, , . . . , (x + i^) ' 

0„(x) = v/F(^W^ [0,,, . . . , 0„](x), (2.98) 

Wc[(;/>di, ■ ■ ■ ,0dj(a:) 
Wc[0di,...,0dj(a; + 1) 
Wc[(/)d,,...,(/)d,,(/>^](a: + l) 
Wc[0.„...,0.„0;.](x) ' ^ • ^ 

■Wc[(^di,...,0dj(x + 1) 



rdQM: B{x) = ^ B{x + l)D{x + I + I) 



D{x) = ^B{x-l)D{x^ 



^ Wc[0d,,...,0rf,,0^](x- 1) 
Wc[(/)di, • • • ,</>d,,</>/.](2;) 



^ -^^^^ Jb{x) - e\/D{x), A^ = JB{x) - JD{x) 6"^ (2.101) 



, def ^Ui=iBix + k-l)Dix + k) 

p [xj — 



Wc [(^di , . . . , 0dj (x) Wci^di , • • • , 0dj (x + 1) ' 
0„(x) = i-lY^/F{xjWc[(l)d,, . . . , 0d„ 0n](2:), (2.102) 
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and 

oQM, dQM: {^n.K) = X{{£{n) - £{dj)) ■ {ct)nAm) {n,meZ>o\V). (2.103) 

It should be emphasised that the Hamihonian as well as the eigenfunction are 
symmetric with respect to di, . . . ,d£, and thus they are independent of the order of {dj}. 
In order to guarantee the positivity of the norm (12.1031) of all the eigenfunctions {0„(a;)} of 
the modified Hamiltonian, the set of deleted energy levels V = {di, . . . , de} must satisfy the 
necessary and sufficient conditions [2S1 [22] 

e 

Y[{m-dj)>0, meZ>o. (2.104) 
i=i 

And the Hamiltonian H is non-singular under this condition. The Crum's theorem in § 12.21 
corresponds to the choice {di, ^2, • • • , di} = {0, 1, . . . , ^ — 1}. 

Starting from an exactly solvable Hamiltonian, one can construct infinitely many vari- 
ants of exactly solvable Hamiltonians and their eigenfunctions by Adler's and Garcia et al's 
methods [25j. The resulting systems are, however, not shape invariant, even if the starting 
system is. For the dQM with real shifts, see the recent work [26J and a related work [55] . 

2.4 Shape Invar iance 

Shape invariance [30] is a sufficient condition for the exact solvability in the Schrodinger 
picture. Combined with Crum's theorem [27], or the factorisation method [IQ] or the so- 
called supersymmetric quantum mechanics [HI [51], the totality of the discrete eigenvalues 
and the corresponding eigenfunctions can be easily obtained. It was shown by the present 
authors that the concept of shape invariance worked equally well in the dQM with the 
pure imaginary shifts [T2[ |2T] as well as the real shifts [20j, providing quantum mechanical 
explanation of the solvability of the Askey scheme of hypergeometric orthogonal polynomials 
in general. Throughout this review, we concentrate on the simple situations that the entire 
spectra are discrete. When continuous spectra exist, however, shape invariance and Crum's 
theorem fail to provide the full set of eigenvalues and eigenfunctions. 

In many cases the Hamiltonian contains some parameter(s), A = (Ai, A2, . . .). Here we 
write the parameter dependence explicitly, 'H(A), ^(A), S{n;X), 0„(x;A), Pn{ri{x; X); X) 
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etc, since it is the central issue. The shape invariance condition with a suitable choice of 
parameters is 

^(A)^(A)^ = kA{X + 5)U(A + 6)+ £{1; A), (2.105) 

where k is a real positive parameter and 6 is the shift of the parameters. In other words 
and ^[^1 have the same shape, only the parameters are shifted by S. The s-th step 
Hamiltonian 7{W in § 12.21 is 7/'*^ = k^'H{\ + s5) + S{s] A). The energy spectrum and the 
excited state wavefunction are determined by the data of the groundstate wavefunction 
(poi^x] A) and the energy of the first excited state £"(1; A) as follows [3ll [121 EDI l2T| : 

n-l 

£(n;A) = ^/s:^^(1;AW), AW = A + s5, (2.106) 

0„(x; A) oc ^(AM)t^(AW)U(A[2])t . . . ^(A["-^')t0o(x; A^). (2.107) 

The above formula for the eigenfunctions 0„(x; A) can be considered as the universal Ro- 
drigues formula for the Askey scheme of hypergeometric polynomials and their g-analogues. 
For the explicit forms of the Rodrigues type formula for each polynomial, one only has to 
substitute the explicit forms of the operator ^(A) and the groundstate wavefunction 4>o{x; A). 
For the nine explicit examples given in fl2.30p - fl2.39p . it is straightforward to verify the shape 
invariance conditions fl2.105p and the energy f l2.106p and the eigenfunction f l2.107p formulas. 

In the case of a finite number of bound states, e.g. the Morse potential, the eigenvalue 
has a maximum at a certain level n, S{n; A). Beyond that level the formula fl2.106p ceases to 
work and the Rodrigues formula f l2.107p does not provide a square integrable eigenfunctions, 
although (^ > continues to satisfy the Schrodinger equation with S{m] A). 

The above shape invariance condition f l2.105p is equivalent to the following conditions: 

oQM: {d^w{x;X)Y -dlw{x;X) = {d^w{x;X + 6)Y + dlw{x;X + 6) + S{l;X), k = 1, 

(2.108) 

idQM: V{x -q;X)V*{x -q;X) = k'^V{x;X + 5)V*{x -ir,X + 5), (2.109) 
V{x + q; A) + V*{x - q; A) = k{V{x; X + 6) + V*{x; X + 6)) - S{1; A), 

(2.110) 

rdQM: B{x + l;X)D{x + l;X) = B{x;X + S)D{x + l;X + S), (2.111) 
B{x; A) + D{x + 1; A) = k{B{x; X + d) + D{x; X + S))+ £{1; A). (2.112) 

For the idQM, the first condition (12.1091) is multiplicative. If Vi, V{, ki and V2, V2, ^2 
satisfy the condition independently, then V = V1V2, V* = V{'V2, k = k\K2 also satisfies it. 
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The second conditions (I2.110p provides a substantial constraint. The situation is the same 
in rdQM. 

It is straightforward to verify the shape invariance for the nine examples fl2.30p - fl2.39p in 
§2. II with the following data: 

oQM: H: A = (null), 6 = (j), k = 1, A) = 2n, (2.113) 

L: A = ^, 6 = 1, K=l, 8{n;\)=An, (2.114) 
J: A=(^,/i), (5 =(1,1), K = l, S{n]X)=An{n + g + h). (2.115) 

It should be stressed that the above shape invariant transformation A — t- A + (5, — )■ 

for L and J, that is, (7— t-^' + I,/;,— t-Zi + I, preserves the monodromy fl2.33p at the regular 

singularities. 

idQM: MP: A = a, 6 = ^, n = 1, £{n;X) = 2n, (2.116) 
W: A = (01,02,03,^4), <5 = (|,|,|,|), K = 1, 

^(n;A) =4n(n + 6i -1), (2.117) 
AW: = (01,02,03,04), (5 = (i, i, 1, i), K = q~^, 

£:(n;A) = (g-"-l)(l-M"-'), (2.118) 

rdQM: M: A=(/3,c), (5 = (1,0), k = 1, £{n;X) = n, (2.119) 

R: X = {a,b,c,d), (5 = (1, 1, 1, 1), k = 1, ^(n; A) = 4n(n + J), (2.120) 

qR: = {a,b,c,d), 5 = (1,1, 1,1), k = 

S{n;X) = {q-^-l){l-dq^~'), (2.121) 

where q(-^^'-^'^'---) = (g-^i , q'^^^ . . .). For more complete lists of exactly solvable oQM see [IHl [51] 
and for idQM and rdQM see [21] and |20] . It is well known that the solutions of Schrodinger 
equations provide those of the corresponding Fokker-Planck equations [56], which describe 
the time evolution of probability density functions. This connection can be readily gener- 
alised to discrete Schrodinger equations, which would correspond to, for example, discretised 
stochastic processes, the Markov chains. The birth and death processes [5l [57] are the best 
known examples. It is interesting to point out that all the exactly solvable rdQM examples 
in [20] and their modifications [261 [H] provide also exactly solvable birth and death processes 

[sg. 

The simplest example of rdQM, the Charlier polynomial, has no shiftable parameter. 
Its shape invariance relation AA^ — A'^A = 1 gives another realisation of the oscillator 
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algebra. Likewise, the simplest example of idQM, the g-Hermite (qH) polynomial has no 
parameter other than q. It is obtained as a special case of the Askey- Wilson (AW) fl2.36p 
by setting aj = 0, j = 1, ... ,4. Like the harmonic oscillator (Hermite) case, it is shape 
invariant and the shape invariance relation {k = q'^) fl2105|) becomes the q-oscillator algebra 
AA^ - g-UU = - 1 itself [19]. 

The shape invariance and the Crum's theorem imply that <pn{x', A) and A + 5) 

are mapped to each other by the operators ^(A) and ^(A)''": 

1 



A{X)(j)n{x; A) = /„(A)0„„i {x; A + <5) x 

1 



^(A)Vn-i(a;; A + 5) = 6„_i(A)0„(x; A) 



X 



: oQM, idQM 
: rdQM 

: oQM, idQM 



(2.122) 
(2.123) 



^B{0;X) :rdQM 

Here the constants /n(A) and 6„_i(A) depend on the normalisation of 0„(x; A) but their 
product does not. It gives the energy eigenvalue. 



Sin; A) = /„(A)6„_i(A). 



(2.124) 



The factor a/ -8(0; A) for rdQM is introduced for later convenience. For our choice of 0n(x; A) 
and Pn{j]{x; A); A), the data for /„(A) and 6„_i(A) are: 




oQM: /„(A) = 
idQM: /„(A) = 
rdQM: /„(A) =f(n;A), 6„_i(A) = 1. 



1 : H 

-2n : LJ ' 



(2.125) 



MP 

W , hn-l{\) 

AW 



n 



MP 

W 

AW 



(2.126) 
(2.127) 



By removing the groundstate contributions, the forward and backward shift operators acting 
on the polynomial eigenf unctions, J^(A) and -B(A), are introduced: 



J-(A) = 0o(a;; A + 5)-^ o ^(A) o 0o(a;; ^) x { ^ 



: oQM, idQM 

v/S(0; A) :rdQM 



A. 

' dr) 



: oQM 
2P) : idQM , 

S(0; A)(/7(x;A)-i(l -e^) : rdQM 



BiX) 



def 



x; o A{Xy o (f)fi{x; X + S) X 



: oQM, idQM 
: rdQM 



(2.128) 
(2.129) 

(2.130) 
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: oQM 



-i{V{x] A)e2P - V*{x; X)e~2P)^{x) : idQM 
^^ {B{x;X)-D{x-X)e'9)ip{x-X) : rdQM 

where c^, ci(?7,A) and 02(^7) are 



(2.131) 




h- g- ig + h + l)r] 




and the auxihary functions if{x) are 



idQM : ip{x) = { 2x 

2smx 



MP 

W 

AW 



^ 1 

2x + d+l 



rdQM : ip{x; X) = < 



d+1 
l-dq 



M 
R 

qK 



(2.133) 



Then the above relations fl2.122p -f l2.123p become 

:F{X)Pn{x- X) = UX)Pn-i{x- X + 5), 

B{X)Pn-l{x- X + 5) = hn-l{X)Pn{x- A), 

where we have used the notation 

Pn{x- X) = PMx;Xy,X). 



(2.134) 
(2.135) 

(2.136) 



Corresponding to (12. 5p . the forward and backward shift operators give a factorisation of the 
similarity transformed Hamiltonian f l2.24p . 



niX) = i3(A)J^(A). 



2.5 Solvability in the Heisenberg Picture 



(2.137) 



As is well known the Heisenberg operator formulation is central to quantum field theory. The 
creation/annihilation operators of the harmonic oscillators are the cornerstones of modern 
quantum physics. However, until recently, it had been generally conceived that the Heisen- 
berg operator solutions are intractable. Here we show that most of the shape invariant dQM 
Hamiltonian systems are exactly solvable in the Heisenberg picture, too [HI [TB] . To be more 
precise, the Heisenberg operator of the sinusoidal coordinate rjlx) 



e**^r?(x)e' 



(2.138) 



21 



can be evaluated in a closed form. It is well known that any orthogonal polynomials satisfy 
the three term recurrence relations [3l [5] 

r^Pn{v) = Ar,Pn+M + BMv) + CnPn^M (n>0), (2.139) 

with P_i{jf) = 0. Here the coefficients A^, Bn and C„ are real and A^-iCn > (n > 1). 
Conversely all the polynomials starting with degree and satisfy the above three term 
recurrence relations are orthogonal (Favard's theorem [59]). These relations can also be 
considered as an eigenvalue equation, in which rj is the eigenvalue. For the polynomials in 
rdQM with the universal normalisation P„(0) = 1 fl2.22p . the coefficients of the three term 
recurrence relations are restricted by the condition 

rdQM: Bn = -{An + Cn) (n = 0,l,...). (2.140) 

For the factorised quantum mechanical eigenfunctions f l2.20p . these relations mean 

T]{x)(j)n{x) = An^n+l{x) + B^M^) + C„0„_i(x) {u > 0). (2.141) 

In other words, the operator ri{x) acts like a creation operator which sends the eigenstate 
n to n + 1 as well as like an annihilation operator, which maps an eigenstate r?, to n — 1. 
This fact combined with the well known result that the annihilation/creation operators of 
the harmonic oscillator are the positive/negative frequency part of the Heisenberg operator 
solution for the coordinate x is the starting point of this subsection. As will be shown 
below the sinusoidal coordinate rj{x) undergoes sinusoidal motion f l2.144p . whose frequencies 
depend on the energy. Thus it is not harmonic in general. To the best of our knowledge, the 
sinusoidal coordinate was first introduced in a rather broad sense for general (not necessarily 
solvable) potentials as a useful means for coherent state research by Nieto and Simmons [50] • 
The sufficient condition for the closed form expression of the Heisenberg operator f l2.138p 
is the closure relation 

[H, [n, vix)] ] = r]{x) Roin) + [n, v{x)] Ri{W) + R-i{Wj. (2.142) 

Here the coefficients Ri{y) are polynomials in y. It is easy to see that the cubic commu- 
tator \H,\H,\H,ri{x)\\\ = (ad'HYri^x) is reduced to ij^x) and [Hjrji^x)] with Ti depending 
coefficients: 
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7]{x) Ro{n)Ri{n) + [n, vix)] [Riinf + i?oCH)) + /?_iCH)i?iCH), (2.143) 



in which the definition {a.d'H)X =^ [HjX] is used. It is trivial to see that all the higher 
commutators (ad "H)"" 77(0;) can also be reduced to ri{x) and \H,ri{x)] with % depending coef- 
ficients. The second order closure fl2.142p simply reflects the Schrodinger equation, which is 
a second order differential or difference equation. Thus we arrive at 

e**«r/(x)e-^*^ = ^(adH)Xx) 

n=0 

pia^{'H)t pia^{'H.)t 

= [^'^(^)] — rm ruv - R'lWM-H)-' 

This simply means that r}{x) oscillates sinusoidally with two energy- dependent "frequencies" 
a±{l-L) given by 

«±(H) = I (i?i(H) ± ^/Ri{V-Y + 4i?o(H) ) , (2.145) 
a+{n) + a^{H) = Ri{n), a+{H)a.{H) = -Rq{H). (2.146) 

The energy spectrum is determined by the over-determined recursion relations £{n + 1) = 
£{n) + a^{£{n)) and £ {n — 1) = £ {n) + a^{£{n)) with £^(0) = 0. It should be stressed that 
for the known spectra {£{n)} determined by the shape invariance, the quantity inside the 
square root in the deflnition of a±{'H) f l2.145p for each n: 

Rr{£{n)f + AR^{£{n)) 

becomes a complete square and the the above two conditions are consistent. For oQM, 
the Hamiltonian and the sinusoidal coordinate satisfying the closure relation f l2.142p are 
classifled and then the eigenf unctions have the factorised form f l2.20p |14] . For dQM we 
assume fl2.20p . The annihilation and creation operators a*^^^ are extracted from this exact 
Heisenberg operator solution: 

e**^r7(x)e-'*^ = a(+)e*"+(^)* + a^-^e'"-^^)* - i?_i(H)i?o(^)"\ (2.147) 
^(±) ±(^[^^^(3,)] _ + i?_i(H)i?o(H)-^)a^(H)) (a+(H) - a_(H))"' 

= ±(a+(?/) -a_(?/))~'([-H,77(x)] +«±(?^)(?7(x) + /2_i(7/)i?oCH)-')), (2.148) 
a(+)t = a(-\ a(+Vn(a;) = A„(^„+i(a;), a^-^ (l)n{x) = Cn<pn-i{x). (2.149) 
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It should be stressed that the annihilation operator and creation operators are hermitian 
conjugate of each other and they act on the eigenstate (12. 1491) . Simple commutation relations 

[•H,a(±)] = a(±)a±(?{), (2.150) 

follow from (12.1481) and (I2.142p . Commutation relations of a^^^ are expressed in terms of the 
coefficients of the three term recurrence relation by (I2.149P : 

^ [a(-),aW]0„ = (A„a+i-^n-iC„)0n. (2.151) 
These relations simply mean the operator relations 

aWflW = /(-H), a^a^-) = g{H), (2.152) 

in which / and g are analytic functions of Ti explicitly given for each example. In other 
words, "H and a*^^-* form a so-called quasi-linear algebra [61]. Various dynamical symmetry 
algebras associated with exactly solvable QM, including the g-oscillator algebra, were explic- 
itly identified in [20| |2T]. It should be stressed that the situation is quite different from those 
of the wide variety of proposed annihilation/creation operators for various quantum systems 
[62], most of which were introduced within the framework of 'algebraic theory of coherent 
states', without exact solvability. In all such cases there is no guarantee for symmetry rela- 
tions like (I2.152p . The explicit form of the annihilation operator (I2.148P allows us to define 
the coherent state as its eigenvector, a^'^tp^ajx) = a'ip{a,x), a G C. See [HI |2T] for various 
coherent states. 

The excited state wavefunctions {(/)„(z)} are obtained by the successive action of the 
creation operator a^~^^ on the groundstate wavefunction (f)o{x). This is the exact solvability 
in the Heisenberg picture. 

The data for the three examples in oQM (I2.30p - (l2.32p are: 

H: i?i(2/) = 0, i?o(?/) = 4, R^i{y)=0, (2.153) 
L: Ri{y) = 0, Ro{y) = lQ, R^i{y) = -8{y + 2g + 1), (2.154) 
J: Ri{y) = 8, Ro{y) = lQ{y+{g + h)^-l), R.,{y) = 16{g - h){g + h - 1). (2.155) 

The data for the three examples in idQM (I2.34p - (l2.36p are: 

MP: Riiy)=0, Roiy)=^ i?_i(y) = 0, (2.156) 
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W: Ri{y) = 2, Ro{y)=4y + hih-2), R_,{y) = -2y^ + (bi - 2b2)y + {2 - bi)h, 
b2 = ^ a^afc, 63 = ^ ajttkai, (2.157) 

AW: R,{y) = {q-h - qhfy\ R,(^y) = [q-\ - q\f [y'^ - [I + q-^fb,) ^ 

R^M = - qh'iih + q~%)y - (l + q''){h + q'%K)) - Ro{y), 

4 

y' = y + 1 + q~^b4, bi = ^aj, 63= ^ ajakai. (2.158) 

j=l l<j<k<l<4: 

The closure relation and the creation/annihilation operators of the MP polynomials were 
obtained in 2001 in |47j. The data for the three examples in rdQM f l2.37p - fl2.39p are: 

M: RM=0, Roiy) = l, R^iiy) = -]^y - Y^^^ (2-159) 
R: Ri{y)=2, Ro{y) = 4y + d^-l, 

R-iiy) = 2y'^ + {2{ab + bc + ca) - {1 + d){l + d))y + abc{d - 1), (2.160) 
qR: R,{y) = {q~\ - qhfy\ R,^y) = ^q'\ - qhf {jj'^ ^ ^q~h + qhfd) , y'"^ y + I + d, 

R^i{y) = {q-^ - q^ f{{l + d)y''^ - [a + b + c + d + d + {ab + be + ca)q'^)y' 

+ (l-a)(l-6)(l-c)(l-dV) 

+ {a + b + c-l-dd+{ab + bc + ca)q-^){l + d)Y (2.161) 

For oQM, the necessary and sufficient condition for the existence of the sinusoidal coor- 
dinate satisfying the closure relation fl2.142p is analysed in Appendix A of [14J. It was shown 
that such systems constitute a sub-group of the shape invariant oQM. We also mention that 
exact Heisenberg operator solutions for independent sinusoidal coordinates as many as the 
degree of freedom were derived for the Calogero systems based on any root system |16j . 
These are novel examples of infinitely many multi-particle Heisenberg operator solutions. 

There were attempts to relate g-oscillator algebras to the difference equation of the q- 
Hermite polynomials [53] • None of them is based on a Hamiltonian, thus hermiticity is 
not manifest and the logic for factorization is unclear. There are several derivations of 
annihilation/creation operators together with the dynamical symmetry algebras (including 
the g-oscillator) acting on the poljTiomials. 
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2.6 Dual Polynomials in rdQM 

The dual polynomials are an important concept in the theory of orthogonal polynomials of 
a discrete variable [2], that is the orthogonal polynomials appearing in rdQM. We show that 
the dual polynomials arise naturally as the solutions of the original eigenvalue problem, the 
Schrodinger equation fl2.17p or f l2.23p . obtained in an alternative way [20]. We do believe this 
derivation of the duality is more intuitive than the existing ones [321 133]. The Hamiltonian 
"H fl2.12p for rdQM is a real symmetric tri-diagonal (Jacobi) matrix of a finite {x^^x = N) 
or infinite (xmax = oo) dimension: 

'^x,y — '^y,xi (2.162) 



n^^y = -^/B{x)D{x + l)6,+,,y - ^/B{x-l)D{x)6,^i,y + {B{x) + D{x))6,^y. (2.163) 



It is well-known that the spectrum of a Jacobi matrix is simple (12. 4p . that is no degener- 
acy. The factor A"^ (I2.10p is a lower triangular matrix with the diagonal and sub-diagonal 
entries only, and A is upper triangular having the diagonal and super-diagonal entries only. 
Throughout this paper we adopt the (standard) convention that the (0, 0) element of the 
matrix is at the upper left corner. The similarity transformed Hamiltonian "H (I2.27P is again 
tri-diagonal: 



^ — ('^x',?/)o<a:,?y<a;max5 '^x,y — B{x){Sx,y — S^+l^y) + D{x){5x,y — (^x-l,j/)- (2.164) 



The similarity transformed eigenvalue problem 'Hv{x) = Sv{x) (I2.23P can be rewritten 
into an explicit matrix form with the change of the notation v{x) — )■ ^{Qo, Qi, . . . , Qx, ■ ■ •) 



Because of the tri-diagonality of "H, the above eigenvalue equations are in fact the three term 
recurrence relations for {Qx} as polynomials in S: 



SQxiS) = B{x){Qxi£)-Qx+ii£))+Dix){Qxi£)-Qx^ii£)) (x = 0, 1, . . . ,x^ax). (2.166) 



Starting with the boundary (initial) condition Qq = 1, Qx{£) is determined as a degree x 
polynomial in It is easy to see 




(2.165) 



g,(0) = l (x = 0,l,.... 



(2.167) 
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When E is replaced by the actual value of the n-th eigenvalue £-{n) (I2.28P in we 
obtain the explicit form of the eigenvector 

). (2.168) 

y=o 

In the finite dimensional case, {Qo{£), . . . ,Qn{S)} are determined by (12.1661) for x = 
0, . . . ,N — 1. The last equation 

SQjviS) = D{N){Q^{S) - Qn-i{S)), (2.169) 

is the degree + 1 algebraic equation (characteristic equation) for the determination of all 
the eigenvalues {£{n)}. 

We now have two expressions (polynomials) for the eigenvectors of the problem (12.231) 
belonging to the eigenvalue £{n); Pn{r]{x)) and Qx{£{n)). Due to the simplicity of the 
spectrum of the Jacobi matrix, they must be equal up to a multiplicative factor a„, 

Pn{r]{x)) = anQx{£{n)) (x = 0, 1, . . . ,Xmax), 

which turns out to be unity because of the boundary (initial) condition at x = (I2.2ip . 
(E22]), (121671) (or at n = (EUl), (E^Z]), (EW])); (n^ax = w) 

Pn{r]m = Pn(0) = 1 = Q^{S{n)), (n = 0, 1, . . . , n^^x), (2.170) 

Po{r]{x)) = l = QM = Qx{£m. (x = 0,l,...,x^ax). (2.171) 

We have established that two polynomials, {Pn{v)} ^"^^ its dual polynomial {Qxi^)}-, 
coincide at the integer lattice points: 

Pn{ri{x)) = Q^{£{n)) (n = 0, 1, . . . , nmax ; = 0, 1, . . . , Xmax)- (2.172) 

The completeness relation, which is dual to the orthogonality relation fl2.50p . 

^max ^max -t 

J2 dlPM^))PMy)) = E dlQx{S{n))QyiSin)) = — (2.173) 

n=0 n=0 ^ 

is now understood as the orthogonality relation of the dual polynomials Qx{S), and the 
previous normalisation constant is now the orthogonality measure. 

The real symmetric (hermitian) matrix T-Lx,y fl2.163p can be expressed in terms of the 
complete set of the eigenvalues and the corresponding normalised eigenvectors 

''T'max 

'^x,y = Yl ^i^)M^)My), (2.174) 

n=0 
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(I)n{x) = dn(f)o{x)Pn{T]{x)) = dn(t)Q{x)Q {jl)) . 



(2.175) 



The very fact that it is tri-diagonal can be easily verified by using the difference equation 
for the polynomial Pn{'r]{x)) or the three term recurrence relations for Qx{£{n)). 
Here is a list of the dual correspondence: 



The functions B{x) and D{x) govern the difference equation for the polynomials Pn{v)i 
the solution of which requires the knowledge of the sinusoidal coordinate rj{x). The same 
quantities B{x) and D{x) specify the three term recurrence of the dual polynomials Qx{£) 
without the knowledge of the spectrum £{n). It is required for them to be the eigenvectors of 
the eigenvalue problem (12.231) . Likewise, and C„ in (12.1391) specify the polynomials Pnij]) 
without the knowledge of the sinusoidal coordinate. As for the dual polynomial Q^(£(n)), An 
and Cn provide the difference equation (in ra), the solution of which needs the explicit form 
of £{n). Let us stress that it is the eigenvalue problem (I2.23P with the specific Hamiltonian 
(I2.12P that determines the polynomials Pn{v) ^-^id their dual Qx{£), the spectrum S{n), the 
sinusoidal coordinate ri{x) and the orthogonality measures (poix)"^ and d?^. 

By comparing the explicit Heisenberg operator solution for r]{x) with the three term 
recurrence relations (12.1391) . the coefficients An and C„ are determined: 



X -H- n, r]{x)^£{n), t]{0) = ^ £{0) = 0, 

B{x)^-An, D{x)i^-Cn, 




(2.176) 



(2.177) 



A 



-0 



i?_i(0)i?o(0)-\ Co = 0, 



(2.178) 



A 



'n 



■n 



R.mn)) + ^(1) (g(n) - B{Q))a^{£{n)) 
«+(£(n))(«+(£:(n))-a_(£(n))) 

R-i{£{n)) + n{l)[£{n) - E(0))a_(g(n)) 
a.{£{n)){a.{£{n)) - a^{£{n))) 



(2.179) 



(2.180) 



There is an important relation among the four important quantities: 



^f(l) + S(0)r7(l) = 0. 



(2.181) 



For the details of the derivation, see (4.45)-(4.53) of [20] 
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2.7 Dual Closure Relation in dQM 

The dual closure relation has the same form as the closure relation fl2.142p with the roles of 
the Hamiltonian "H and the sinusoidal coordinate ri{x) interchanged [20] : 

[V, [v, n]] = n i?;J"^'(r/) + [ri, n] i?f ^^(r]) + i?^T'(r/), (2.182) 

in which 

i?f ^'(r/(z)) = {r]{x - i(3) - ri{x)) + {r]{x + i(3) - ri[x)) , (2.183) 

i?;J"^'(r/(x)) = -{ri[x-%^)-ri[x)){r]{x^i^)-ri[x)), (2.184) 

i?^T'(r/(x)) = e{y^{x) + V_{x))Rf-''\^{x)). (2.185) 

The dual closure relation is the characteristic feature shared by all the 'Hamiltonians' 
"H which map a polynomial in riix) into another. Therefore its dynamical contents are not 
so constraining as the closure relation, except for the rdQM exactly solvable case, where 
the closure relation and the dual closure relations are on the same footing and they form a 
dynamical symmetry algebra which is sometimes called the Askey- Wilson algebra [511 EH 
[331 [20]. 

2.8 Bochner's Theorem 

In 1884 [65], Routh showed that polynomials satisfying the three term recurrence relations 
and a second order differential equation were one of the classical polynomials, the Hermite, 
Laguerre, Jacobi and Bessel. Later in 1929 [66], Bochner classified all polynomial solutions 
to second order Sturm-Liouville operators with polynomial coefficients and arrived at the 
same conclusions. See §20.1 of [[5] for more details. This was a kind of No-Go theorem in 
oQM, since it declared that no essentially new exactly solvable oQM could be achieved as 
the solutions of the ordinary Schrodinger equation. Thus avoiding Bochner's theorem was 
one of the strongest motivations for the introduction of the discrete quantum mechanics, 
whose difference Schrodinger equations were not constrained by the theorem and provided 
various exactly solvable examples [201 EI] • As shown in the new orthogonal (exceptional 
(Xi) Laguerre and Jacobi) polynomials [23, 15S1 IS] were discovered in an attempt to evade 
the restrictions of the theorem by allowing the polynomials to start at degree i > 1. The 
new orthogonal polynomials in dQM, the exceptional (X^) Wilson, Askey- Wilson, Racah 
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and g-Racah polynomials satisfying second order difference equations, were constructed by 
the present authors in less than two years. 

Here we present Bochner's theorem for dQM. As in the original Bochner's theorem, 
the conditions for the lowest three degrees n = 0, 1 and 2 are essential to determine the 
constraints. Reformulation of Bochner's theorem for second order difference equations was 
pursued by several authors. In |69], it was shown (in our language) that if orthogonal poly- 
nomials in = cosx satisfy q'-difference equations f l2.23p . then the conditions characterise 
the Askey- Wilson polynomials. In [70] , it was shown that if the g-difference equations fl2.23p 
have polynomial solutions, then the sinusoidal coordinate r]{x) is at most g-quadratic and 
that the polynomials are at most Askey- Wilson polynomials. In these papers, the distinction 
between the pure imaginary and the real shifts is blurred. 

The starting point is the difference equation for the polynomials fl2.23p with fl2.26p : 

eV4x) {Pnivix - - Pnivix))) + eV.{x) {Pn{v{x + t/3)) - Pnivix))) 

= £in)PMx)) (n = 0,l,...), (2.186) 

and the three term recurrence relations fl2.139p . The above equation is trivially satisfied for 
n = 0, since S{0) = 0. Two relations for n = 1 and n = 2 are linear equations 

PM^ - 2/3)) - Pi(?7(x)) P^{r]{x + t/3)) - PMx))\ (V+{x)\ ^ ( e-^£{l)PMx)) 
Pnivix - m - P2ivix)) P2Hx + ip))-P2{ri{x))) \y^{x)) \e-^£{2)P2{r]{x)) 

which determine the potential functions V+{x) and V-{x) uniquely. By using the explicit 
forms of Pi and P2, 

AoPi(?7) = - 5o, AoAiP2(?7) = (r/ - So)(?7 - Pi) - ACi, (2.187) 

derived from the three term recurrence relations (12.1391) . the above linear equation gives 

y^^^) ^ S2 + Si7]{x±iP) ^ ^2.188) 

{ri{x =1= — r]{x)) {rj{x =F "iP) — v{x ± ^/3)) ' 

in which 6*1 and 5*2 are a linear and a quadratic polynomial in r]{x), respectively: 

S, = -s-'£il)AoPMx)), 

S2 = e-'S{2)AoA,P2{v{x)) - e-'S{l)AoPMx)){AoPMx)) - Pi). (2.189) 

These V± have essentially the same forms as those which will be presented in §[3l If the 
symmetric shift-addition property of the sinusoidal coordinate rji^x) (13. 4p is satisfied, the 
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above expressions for V±{x) fl2.188p are equivalent to V± in fl3.2l) - fl3.3p with L = 2. Moreover 
when the symmetric shift-muhiphcation property fl3.5p is satisfied, these potential functions 
give the exactly solvable models. For rdQM, the sinusoidal coordinates satisfying (13.40 can 
be classified into five types fl3.15l) - fl3.19l) [20] and they also satisfy (13. 5p . 

3 Unified Theory of Exactly Solvable dQM 

In this section we present a simple theory of constructing exactly solvable 'Hamiltonians' 
in dQM based on the two required properties of the sinusoidal coordinates (I3.4p -( P3]) . The 
general strategy is to construct the similarity transformed 'Hamiltonian' "H (I2.26P in such a 
way that it maps a polynomial in r]{x) into another: 

?^V„ C V„+L_2 C Voo (nG Z>o), (3.1) 

where L is a fixed positive integer. Here V„ is defined by (I2.29P and Voo == lini„^oo V„. 
When L = 2, the above relation (13. ip is simply the lower triangularity of the 'Hamiltonian' 
H, leading to exact solvability. 

In the following we will take the similarity transformed Hamiltonian "H (I2.24p instead 
of "H as the starting point. That is, we reverse the argument and construct directly the 
'Hamiltonian' Ti (12.261) based on the sinusoidal coordinate r]{x). This section is a brief 
review of [23j . 

3.1 Potential Functions 

The general form of the 'Hamiltonian' "H mapping a polynomial in ri{x) into another is 
achieved by the following form of the potential functions V±(x): 

V^{x) = , (3.2) 

V±{x) = r/(x)S(x T ^/3)^ (3.3) 

k,l>0 
k+l<L 

where L is a natural number indicating the degree of rj{x) in V±{x) and Vk,i are real constants, 
with the constraint ^ f| ^ 7^ 0. It is important that the same Vk^i appears in both V±{x). 

k+l=L 

The 'Hamiltonian' H with the above V±{x) maps a degree n polynomial in ri{x) to a degree 
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n + L — 2 polynomial. This can be shown elementarily based on the two basic properties of 
the sinusoidal coordinates called the symmetric shift- addition property: 

r]{x - i^) + r]{x + ip) = (2 + r^^\{x) + r% (3.4) 

and the symmetric shift-multiplication property: 

r]{x — ij3)rj{x + i(5) = {r]{x) — T]{—if3)) {r]{x) — r](i(3)), (3.5) 

together with r]{x) ^ r^(x — i/3) ^ r]{x + i/?) ri{x). Here r^j^'^ and r*^] are real parameters. 
In fact these parameters also appear in the three functions Rq, Ri and in the closure 
relation (I27[l2|) : 

Ri{y) = r['^y + rf>\ R,^y) = r^^)y^ + ri^^y + r^^\ R-M = r%' + r^'lv + r^^l (3.6) 
Here are the lists of the known sinusoidal coordinates satisfying the above two conditions 



fl3.4p -( l33|) . There are eight sinusoidal coordinates for the idQM: 

(i) : r}{x) = x, — oo <x < oo, 7 = 1, (3.7) 

(ii) : r]{x) = x^, <x < oo, 7 = 1, (3.8) 

(iii) : ri{x) = 1 — cosx, <x < tt, 7 G M^q, (3.9) 

(iv) : ?7(x) = sinx, -f<x<f, 7 G M^o, (3.10) 

(v) : ?7(x) = l-e"^, -oo<x<oo, 7 G M^o, (3-11) 

(vi) : r]{x) = e''-l, -oo<x<oo, 7 G M^o, (3.12) 

(vii) : ri{x) = cosh X — 1, 0<x<oo, 7 G M^o, (3.13) 

(viii) : rj{x) = sinhx, —00 <x < 00, 7 G M^o, (3-14) 

and five sinusoidal coordinates for the rdQM: {0 < q < 1) 

(i) ' : r](x) = x, (3.15) 

(ii) ' : r]{x) = e'x{x + d), e' = | d < -iV ' '^^'^^'^ 

(iii) ': ?7(x) = l-g^ (3.17) 

(iv) ' : r]{x) = q^"" - 1, (3.18) 

(v) ': r^{x) = e\q'^ -dq% e' = { \ tt^N ■ (3-19) 
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As shown in detail in §4C of [20], the above five sinusoidal coordinates for rdQM f l3.15p - 
flS.igp exhaust all the solutions of fl3.4p -f l33|) up to a multiplicative factor. On the other 
hand, those for idQM (i)-(viii) f l3.7p -( l3TT4|) are merely typical examples satisfying all the 
postulates for the sinusoidal coordinate fl3.4p - fl3.5p and the extra one used for the shape 
invariance, (3.50) of [ 23] . 

The essential part of the formula (13. 2p is the denominators. They have the same form as 
the generic formula for the coefficients of the three term recurrence relations of the orthogonal 
poljTiomials, (I2.179P and (I2.180p (see (4.52) and (4.53) in [20]). The translation rules are 
the duaUty correspondence itself, (1277761) - (121771) (see also (3.14)-(3.18) in [20]): 

S{n)^r]{x), -An^V+{x), -Cn^V-{x), 
aj^{£ {n)) ^ rj{x - iP) - r]{x), a_[S{n)) ^ r]{x + - r]{x). (3.20) 

Some of the parameters Vk,i in (13. 3p are redundant. It is sufficient to keep Vk,i with 
/ = 0, 1. The remaining 2L + 1 parameters Vk,i {k + I < L, I = 0,1) are independent, with 
one of which corresponds to the overall normalisation of the Hamiltonian. See § II of [23j for 
the actual proof of the property (13. ip . 

The L = 2 case is exactly solvable. This corresponds to the general hypergeometric 
equations having at most degree two polynomial coefficients. Since the Hamiltonian of 
the polynomial space H is expressed as an upper triangular matrix, its eigenvalues and 
eigenvectors are easily obtained explicitly. For the solutions of a full quantum mechanical 
problem, however, one needs the square- integrable groundstate wavefunction (poi^) (I2.18p . 
which is essential for the existence of the Hamiltonian "H and the verification of its hermiticity. 
These conditions would usually restrict the ranges of the parameters fo,o, • • • , "^2,0- It is easy 
to verify that the explicit examples of the potential functions V{x), V*{x) and B{x), D{x) 
in (I2.34p - (l2.39p are simply reproduced by proper choices of the parameters {vk,i}- It should 
be stressed that the above form of the potential function (I3.2p - (l3.3p provides a unified proof 
of the shape invariance relation (I2.105p . the closure relation (I2.142p and the dual closure 
relation (I2.182p in the % scheme, see ^3J for more details. This is in good contrast with the 
various explicit examples presented in the preceding section. The two solution methods, the 
shape invariance and the closure relation are verified for each example. 
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3.2 Askey- Wilson Algebra 

Here we briefly conunent on the Askey- Wilson algebra, which are generated by the closure 
plus the dual closure relations. By simply expanding the double commutators in the closure 
fl2.142p and the dual closure fl2.182p relations, we obtain two cubic relations generated by 
the two operators "H and r]: 

H^T] - (2 + r^^^)Hr]n + r^n^ - rf\Hr] + r]H) - r^^\ = r^^\n^ + r^H + r% (3.21) 

ifn - (2 + r^^'^)r]H'n + Hrf - r^^lirjU + Uv) + r]{-i^)r]{i 13)1-1 = rf'^rf' + r^^lr] + evq^q. 

(3.22) 

From its structure, the closure relation is at most linear in rj and at most quadratic in "H. 
So the l.h.s. of f l3.2ip has terms containing one factor of rj and the r.h.s, none. It is simply 
R-i{'H). Likewise, the l.h.s. of f l3.22p has terms containing one factor of l-L and the r.h.s, 
none. It is simply R^i^ijj). In f l3.22p . 'r]{—i(3)ri{i(3) is just a real number, not an operator. 
These have the same form as the so-called Askey- Wilson algebra, which has many different 
expressions [511 IMl ED While the Askey- Wilson algebra has no inherent structure, the 
closure relation fl2.142p has the right structure to lead to the Heisenberg operator solution 
for rjlx), whose positive and negative frequency parts are the annihilation-creation operators 
[131 l20l [2T| . It is the Hamiltonian and the annihilation-creation operators that form the 
dynamical symmetry algebra of the system [201121], not the closure or dual-closure relations, 
nor the Askey- Wilson algebra relations. 

3.3 Quasi-Exact Solvability in dQM 

The unified theory is general enough to generate quasi-exactly solvable Hamiltonians in the 
same manner. The quasi-exact solvability means, in contrast to the exact solvability, that 
only a finite number of energy eigenvalues and the corresponding eigenf unctions can be 
obtained exactly. Many examples are known in the oQM [TH [721 [73], but only a few are 
known in dQM in spite of the proposal that the sl{2,R) algebra characterisation of quasi- 
exact solvability could be extended to difference Schrodinger equations [72j. The unified 
theory also incorporates the known examples of quasi-exactly solvable Hamiltonians in dQM 
[T71 [T8l [22] . A new type of quasi-exactly solvable Hamiltonians is constructed. The present 
approach reveals the common structure underlying the exactly and quasi-exactly solvable 
theories, in particular, the important roles played by the sinusoidal coordinates. 
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The higher L > 3 cases are obviously non-solvable. Among them, the tame non-solvability 
of L = 3 and 4 can be made quasi-exactly solvable (QES) by adding suitable compensation 
terms. This is a simple generalisation of the method of Sasaki & Takasaki [74j for multi- 
particle QES in the oQM. For a given positive integer M, let us try to find a QES 'Hamil- 
tonian' "H, or more precisely its modification "H', having an invariant polynomial subspace 

H'Vm C Vm- (3.23) 
For L = 3, "H' is defined by adding one single compensation term of degree one 

n"^^^n-eo{M)7]{x), (3.24) 

and we have achieved the quasi-exact solvability T-i'Vu ^ Va/. Known discrete QES examples 
belong to this class [HI [22]. 

For L = 4 case, "H' is defined by adding a linear and a quadratic in ri{x) compensation 
terms to the Hamiltonian %: 

n^t^U- eo{M)r]{x)^ - ei{M)r]{x), (3.25) 

and one condition between v^^i and ^4,0 is imposed. Then we have T-L'Vm ^ Vm- This type 
of QES theory is new. For explicit forms of eo(M) and ei(M), see p3] . 

4 New Orthogonal Polynomials 

In this section we present a brief review of the recent hot topic, the discovery of infinitely 
many new orthogonal polynomials satisfying second order differential or difference equations. 
They are obtained as the main parts of eigenfunctions of exactly solvable quantum mechanics 
in all the three categories, oQM, idQM and rdQM. They are the exceptional (Xi) Laguerre 
(XL) and Jacobi (XJ) polynomials in oQM, the exceptional Meixner-Pollaczek (XMP), Wil- 
son (XW) and Askey- Wilson (XAW) polynomials in idQM and the exceptional Meixner 
(XM), Racah (XR) and g-Racah (XgR) polynomials in rdQM. The XMP and XM are new 
results. For any positive integer £ these Xi polynomials have degrees i + n, n = 0,1,2,... 
and are orthogonal with each other (14.43 p . (14.501) . (I4.55P with respect to explicitly given pos- 
itive definite weight functions. Thus they do not satisfy the three term recurrence relations 
(I2.139P and therefore the restrictions of Bochner's theorem and its discrete counterparts 
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in §2.81 do not apply. As will be explained later, the quantum mechanical Hamiltonians for 
these exceptional polynomials are deformations of those for the original polynomials in terms 
of a degree £ eigenpolynomial with twisted parameters. The deformation terms satisfy the 
same type of equations as the original polynomials, but with the twisted parameters. The 
parameters are so chosen that the deformed system retain the shape invariance, too. The 
general scheme of the deformation is common for various exceptional polynomials as shown 
in detail below. The original system corresponds to £ = 0. 

It should be stressed that the Jacobi polynomials provide infinitely many global 
solutions of Fuchsian differential equations having 3 + £ regular singularities [9]. 

4.1 Deforming Polynomials 

In order to introduce new orthogonal polynomials, we prepare the deforming polynomials. 
The deforming polynomials ^^{Tf) are eigenpolynomials Pi{rf} with twisted parameters (or 
twisted coordinate) and they are crucial for the construction of new orthogonal polynomials. 
Two identities relating ^£ with A and A + <5, fHA5D - fH:T6D and fOHD -f lCTll . will 

play important roles in the derivation of various results. 

For oQM there are two types of the exceptional Jacobi polynomials, called XJl and 
XJ2. Their confluent limits give the exceptional Laguerre polynomials, XLl and XL2. The 
deforming polynomials are: 




XLl 



XL2 
XJl 
XJ2 



(4.1) 



Explicitly they are 




(4.2) 



The two identities relating ^^(77; A) and ^^(r/; A + <5) are: 



di(A + l8)Wn; A) + d^d^Um A) = d,{\)Um A + 5), 
ci3(A, A + 5) + ^ d^Uri: A + (5) = 4(A + 15, i)Uvi A), 




(4.3) 



(4.4) 



where di{\), (^2(^7) and d3{\,i) are given by 




XLl 

XL2, XJ2 
XJl 




XLl 
XL2 

XJ1/XJ2 



(4.5) 
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ds{\i)={ 1 



9 + e-l :XL1, XJl 

XL2 . (4.6) 

XJ2 



Note that these two identities, (14. 3 p and (14. 4p . imply the differential equation for the de- 
forming polynomial, 

ci?/,A,^j \ Ci(r7,t(A+(^-l)(5)) :XL2, XJl, XJ2 ' ^^'^^ 

^ ' ' \S{i;t{\)) :XL2, XJl, XJ2 ' ^ ' 

which is to be compared with the differential equation for the eigenpolynomial (I2.23p . 

C2ir])d^^P„iv; A) + ciirj, X)d^Pniv. A) = -l£{n; X)Pn{m A)- (4.10) 

Note that £{t, t(A)) = £{(i; t(A+ (£— 1)<5)). The deforming polynomial iiij]] A) has the same 
sign in the orthogonality domain, that is (0, oo) for L and (—1, 1) for J, see (2.39), (2.40) of 
[35] and (3.2) of [38j. 

For idQM, we restrict parameters: 

XMP : ^ even, (4.11) 
XW:ai,a2GM, {03, 04} = {03, 04} (as a set), 0<aj<Reafc (j = 1, 2; A; = 3, 4), (4.12) 
XAW:ai,a2GM, {og, 04} = {03, 04} (as a set), 1 > > |afc| (j = 1, 2; A; = 3, 4). (4.13) 

The deforming polynomials for XMP, XW and XAW are: 
The two identities relating ie.{r]{x)] A) and ^e{i]{x); X + 6) are 



^{v*ix; A + iS)e-^P - v,{x; A + iS)e"P)Uvix); A) = /,,o(A)6(r/(x); A + S), (4.15) 
^{V2ix; X+ii- l)(5)eiP - v;{x; X + {i- l)(5)e-i^)6(r?(x); X + 6) 

= k,o{X)Uv{xy,X), (4.16) 

where vi{x; A), V2{x; A) are the factors of the potential function V{x; A): 

T/^ \^ ^ Mx; X)v2{x;X) 

^^^'^) = -^ V^(x)^(x-.i)' ^^-'^^ 
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i[a + ix) 
Vi{x;X) = \ Ul=ii(^j + 



XMP fi 

XAW ' -'■■^■r^^ 



tx) 



XMP 
XW . 
XAW 
(4.18) 



The constants fi^W and bi^nW are given by 



2a + n 



dcf 



ai + a2 + n 



[1 - aiaiq'- 



XMP 



XW , 6^^„(A) = <( a3 + a4 + n + 2£- 1 



XAW 



XMP 
XW . 
XAW 
(4.19) 

These two identities f l4.15p and f l4.16p imply the difference equation for the deforming poly- 
nomial. 



{V{x- t(A + {^- 1)6)) {e^^ -l) + V* (x; t(A + {I - 1)6)) (e'^^ - l))eK^(^); A) 



£{t,i{\))Ur,{x)-\), 



which should be compared with the difference equation for the eigenpolynomial 



(4.20) 



{V{x- A)(e^P - 1) + V\x- A)(e-^P - l))F„(r](x); A) = £{n- X)Pn{v{x); A). (4.21) 

Note that £{i;t{\)) = £{i;t{X + {i — 1)<5)). For the appropriate parameter ranges, the 
deforming polynomial ^^(?](z);A) has no zero in the rectangular domain Xi < Rex < X2, 
|Ima;| < I7I, which is necessary for the hermiticity of the Hamiltonian. 
For rdQM, the deforming polynomials for XM, XR and XgR are: 



: XM 



6(x; A) = ^i]{x; X+{i- 1)6); X) , F„(a;; A) = Pn{v{x; A); A) , 

^ X) / + f3 + i-l);X + {e- 1)6) 

\ Pe{x;t{X+ii-l)6)), t(A)^=^f(A4-Ai,A4-A2,A3,A4) :XR, XgR^ 

which satisfies the normalisation 



(4.22) 
(4.23) 

(4.24) 



6(0; A) = 1. 

The two identities relating ^(x; A) and ^(x; A + (5) are [TT] : 

, , \, Uix; X + i6)- <(x; A + i6)e^)Ux; A) = koWUx; A + 5), (4.25) 
(p[x; X + 16 + 6) ^ ' 

, , ] Uix; A + (^ - 1)6) - v^{x; A + (£ - l)6)e-^)Ux; X + 6) 

ip{x] X+ [i - 1)6 + 6) ^ ^ 

= 6,,o(A)6(x;A). (4.26) 
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Here ff (x; A), ^^^(x; A), ff (x; A), vl^i^x] A) are the factors of the potential functions B{x] A) 
and D{x; A): 



B{x;X) 
D{x; A) 



(p{x; A + d)ip{x + |; A + ^) ' 
A)t>^(a;; A) 
A + 6)ip{x — |; A + 5) ' 



— C2 



: XM 
: XR 



1-d 

c^^B{x; A) 
t;f(x;A)'=^J d-\x + c){x + d) 



(1 -ag^)(l-6g^') : XgR 



XM 
XR 



1 - d 

_ 1 

— C 2 



(1 -cg^')(l-(ig^) :XgR 



;f (x- A) =^ < ^ ^x + (i - a)(x + d - b) 



XM 
XR 



A) 



1 - d 
c^D{x] A) 



abd-\l - a-Mg^)(l - : XgR 



XM 
XR 



1-d 



c(l -c"Vg^)(l -g^) :XgR 



where 6 is 

^def f (-1,0) :XM 

\ (0,0,-1,-1) :XR,XgR ■ 

The constants fi^nW and 6£,„(A) are given by 

( l-cP + n + 2i-l 



An (A) =^ <( 



XM 



/3 + ^-l 

(a + 5_^ + ^) ^ + ^^ + ^-^ :XR 
c + i-1 

g-"(l - aM-^g") ^ ^ : XgR 



^^,n(A) 



dcf I 1 - c 



(/3 + £-l) 



1 - cg^-i 



XM 



XR • 

XgR 



(4.27) 
(4.28) 

(4.29) 
(4.30) 
(4.31) 
(4.32) 



(4.33) 



(4.34) 



(4.35) 



c + i-1 
1 - cq'-' 

Again these two identities f l4.25p and fl4.26p imply the difference equation for the deforming 
polynomial, 

{B{x; A)(l - e^) + 5(x; A)(l - e-'))Ux; A) = £{£; X)Ux; A), (4.36) 
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XM 

XR, XgR ' 
XM 

XR, XgR ' 



(4.37) 
(4.38) 



f5.,. 4^f / ~D{-{x + + 1); A + (^ - 1)6) 
^ ''^^ \ B{x;t{X+{e-l)6)) 

~ dcf / -B{-ix + + 1); A + (^ - 1)6) 
^{x,A) - I /)(^.t(A + (£-l)5)) 

£{i;X) = I ^(^I^J;^]) :XR, XgR ' ^^'^^^ 
to be compared with the difference equation for the eigenpolynomial 

{B{x; A)(l - e^) + D{x- A)(l - e-^))P,,{x; A) = ^(n; A)P„(x; A). (4.40) 

Note that £{i;t{X)) = £{i;t{X + {i — 1)6)). For the appropriate parameter ranges, the 
deforming polynomial ie{x] X) is positive at integer points x = 0, 1, ... , x^^^ + 1 (x^^^ =^ 

■^max ^) • 

Note that the range of parameters can be enlarged. For XJ, see [75] . 
4.2 New Orthogonal Polynomials 

In this subsection we present the explicit forms of the exceptional orthogonal polynomials 
and the orthogonality relations. They are bilinear in the deforming polynomial ^£ (with A 
and A + 5) and the eigenpolynomial P„ (with X + £6 + 6). 
For oQM, the exceptional polynomials are: 

PUv, A) = y4TT (^2(^)6(r/; X)dr,Pn{v; X + e6 + 6) 
Je,n{X) ^ 

- d^iX)Ci{v■, X + 6)Pn{v; A + M + 5)) X I • /xj2^ , (4.41) 

where 6 is 

S I • XL1/XL2 

^ ~ \ t(1,-1) : XJ1/XJ2 ' ^ ' 

and fi,n{X) will be given in f l4.45p . They satisfy the orthogonality relation: 

/ ijgix; A)^ Pe^nivix); X)Pe^rn{'n{x); X)dx = he,n{X)Snm, (4.43) 

J X\ 

where the weight function '^^(x; A)^ will be given in fl4.63p . The normalisation constant 
hin{.X) is related to /i„(A) f l2.5ip as 

h,n{X) = h^KiX + 16 + 6)= ^5n(A)/o,n(A + M) 

' An(A) ^ /,,„(A) V(A + M) ^' ^ ^ 

40 



where fe,n{X) and &£,n(A) are given by 



fe,nW = 2x 



-1 


: XLl 


n + g + l 


: XL2 




: XJl 


n + g + l 


: XJ2 



bi^nW = 2 X 



-{n + g + 2i-^) :XL1 

dcf „ 11 



[n + g + 2i-l 



n + h + 2l-\ 



2 



XL2 

XJl ■ 

XJ2 
(4.45) 



The deforming polynomial ^^(?7;A) has the same sign in the orthogonality domain, that is 
(0, oo) for XL and (-1, 1) for XJ. 

The two types of the XJ polynomials are the mirror images of each other 

Pf:^\^; 9. h) = {-lY^-Pll\-x- h, g), (4.46) 

due to the parity property of the Jacobi polynomial Pn°''^\—x) = (— l)"Pi^'"''(x), see (50)- 
(54) of [8J . The Laguerre polynomials are known to be obtained from the Jacobi polynomials 
in a certain limit: 

lim P^"' (1 - 2x/3-^) = L^^^ (±x) . (4.47) 

This connects the XJ polynomials to the XL polynomials as shown in (43), (58) of [S]. 
For idQM, the exceptional polynomials are 



PeAvi^y, A) = . ' ^ U{x; A + id)Uvix + ^i); A)P„(?7(x - 2^); A + M + ^) 
fe,nW(p{x) ^ 

- vlix- A + id)Uvix - ^i); A)P„(?7(x + ); A + M + ^)) , (4.48) 

where d is 

A / i • ™^ (A AQ^ 

^=[ hi-h-l) :XW,XAW • (4-49) 

They satisfy the orthogonality relation: 



X2 



ijje{x; A) P£,„(77(x); A)P£,m(77(x); \)dx = V„(A)5„m, (4.50) 



where the weight function ■0^(x;A)^ will be given in (I4.67p . The normalisation constant 
hi^nW is related to /i„(A) f l2.53p as 

heniX) = h^K{X + l5+~5) = ^5n(A)/o,n(A + M) 

' fi,nW ^ ' /,,„(A)6o.n(A + M) ' ^ ' 

We have verified by numerical calculation for the parameter ranges in (14. 11 1) - (14.131) that 
the deforming polynomial ^^{rj{x)] A) has the same sign in the orthogonality domain, that is 
X e (-00, oo) for XM, x G (0, oo) for XW and x E (0, vr) for XAW. 
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For rdQM, the exceptional polynomials are jllj : 



P,,„(x;A)^=^'F,,„(r/(x;A + M);A) 

dcf 1 1 ( B 



<(x;A + M)6(x;A)P„(x + l;A + M + 5) 

/An(A) ^{X]\ + 16 + 6) 

- v^{x] A + + 1; \)Pn{x- A + M + , (4.52) 

which satisfy the normalisation 

P,,„(0;A) = 1. (4.53) 

Here 6 is 

^defr(-l,0) :XM 

\ (0,0,-1,-1) :XR, XgR • ^ > 

The orthogonahty relation reads (x^^^ = x^ax - ^, ^^Lx =^ 3;max - ^): 

E T^rrf ^)^^-^^' ^) = (n, m = 0, 1, . . . , nLx), (4-55) 

^ 6(1; A) rf^,n(A)^ 

where the weight function '^^(x; A)^ will be given in fl4.7ip . The normalisation constant 
di^nW^ is related to (i„(A)^ fl2.56p as 

d,..ixf ^ d.ix+is+~sr 1^4^ ^ .„(A+M)^ tlfr^^T'^^i^'^ ^ ^'-''^ 

where se{X) is 

f 1-c 1 

: XM 

c f3 + i-l 

(d — a)(d — b) 

^ ' :XR . (4.57) 



ic + £-l){d + i) 

,{l-ar^d){l-_^ 

""^^ ^ (1 - cq^~^)il - dq^) -^^^ 
We have verified by numerical calculation for the parameter ranges in fl2.37p - fl2.39p that 
the deforming polynomial ^(x; A) has the same sign in the orthogonality domain, that is 

T = n 1 

w, X, . . . , ^uiax' 

The exceptional polynomial -P^,n(^; A) is a degree i + n polynomial in r] but has only n 
zeros in (0, oo) for XL, (—1,1) for XJ, {ri{xi),ri{x2)) for idQM and (0,?7(xmax)) for rdQM. 
For £ = 0, the expressions of Pe^niv'i A) reduce to those of the original polynomials. 

These new orthogonal polynomials do not satisfy the ordinary three term recurrence 
relations fl2.139p . For the XL and XJ polynomials, a new type of bi-spectrality is demonstrated 
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in §4 of [36j- The bi-spectrality of the other new orthogonal polynomials, XW, XAW, XR 
and XgR are to be explored. A different type of three term recurrence relations for the XL 
and XJ polynomials is reported in §10 of [9J. The generating functions Yl'^=o t"'Pe,n{V'i ^) 
the XL and XJ polynomials are reported in §9 of [9J. For the XLl and XL2 polynomials, 
the closed form expression of the double generating functions 



oo oo 

def 



G{s,t,ri;\) = }^}_^sTPeAv;^) 

e=o n=0 

are given there, too. The zeros of orthogonal polynomials have always attracted the interest 
of researchers. In [7B] the behaviours of the zeros of the XLl, XL2 and XJ2 polynomials 
Pe,n{v'^ ^) explored as the parameters change. The structure of the extra zeros of the 
other new orthogonal polynomials, XW, XAW, XR and XgR are also an interesting problem. 

4.3 Deformed Hamiltonians 

The new orthogonal polynomials are discovered as the main part of the eigenfunctions of 
quantum mechanical systems in which the original Hamiltonians are deformed in terms of 
the deforming polynomial ^e. The parameters are so chosen that the deformed systems retain 
shape invariance [3 [TOl |35l [HI [371 |36l [11] . 

The deformed Hamiltonian and their eigenfunctions have the following forms: 

niiX)=Ae{XyM\), (4.58) 
ne{X)cPi,nix; A) = ^,,„(A)0,,„(x; A), £,,„(A) = ^(n; A + M), (4.59) 
0,,„(a;; A) = A)P^,n(r/(x; A + id); A). (4.60) 

The deformation is achieved additively at the level of the prepotential for oQM and multi- 
plicatively at the level of the potential functions V{x), V*{x), B[x) and D{x) for dQM. For 
oQM, operators ^^(A), ^^(A)^ and i/jei^x; X) are [7l[8]: 

MX) = -^-d,We{x;X), A(A)t = -A_a,«;,(x;A), (4.61) 

we{x; A) w{x; X + iS) + log ^^(^(^)- ^ + , (4.62) 

^e{7]{x);X) 

def Mx;X + e6) ^ ^^^(,^^) ^ XMvix); X + S). (4.63) 

U{v{x);X} 



For idQM, they are [TO]: 

MX) "=^^(eiyy;(x;A) - e'^V^^l^; A) ) , 
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At{X)^ = -i{^V,{x-\)e'^P - ^Jv*{^)e-^, (4.64) 
Vi{x] A = V[x] X + £6) — 4-— — — , 4.65 

14 [x; X) = V [x; X + id) / c ( ( \ \ ^ a\ ' ^'^^ 

Mx; X) M^-A + ^S) 

V6(^(a; + ^i); A)^^(77(x - i^); A) 



For rdQM, they are [IT]: 



^,(A) =V5^(a:;A)-eVA(x;A), ^,(A)^ = ^^.(x; A) - A) e'^ (4.68) 

S,(^; A) 5(x; A + 6(x + l;A + <5) 

^ ^ ^ ^a^+l;A) e.(^;A + 5) ^ ^ 

D,(x; A) D{x- X + ^s) ^^i^ + ^) - I; X + 5) ^ 

Ci{x; A) ^^(x; A + <5) 

Mx;X) = Mx;X + £S)J-. ^^^2^ . (4.71) 

These potential functions satisfy the boundary conditions ^'^(O; A) = and Be{N — i; A) = 0, 
and the functions (pi^ix; A) and ipii^x; A) satisfy the normahsation conditions (;/>^,o(0; A) = 1 
and '0^(0; A) = 1. 

These deformed systems are shape invariant: 

Ai{X)MXy = kMX + SYMX + S) + ^,,i(A), (4.72) 

or equivalently fl2.108p -( l2.112p with replacements {w,V, B, D) — )• (wi, Vi, Bi, Di). The shape 
invariance of these deformed systems is shown analytically by making use of the two identities 
of the deforming polynomials ^-^^^ ^H^ - ^H^ and (g^SD-dlSSD [2S1EZIIII]. 
The action of the operators A(,{X) and .4^ (A)'^ on the eigenf unctions is 

f 1 : oQM, idQM 

Ae{X)(t)i^n{x\X) = f^^n{X)(t)e,n-i{x]X + 5) ^ I 1 ■ rdOM ' (^-^S) 

A,{X)^,,^.,[x- X + 5)= h,n-i{X)<PeAx; A) x | ^^/^^ [ ^^c^'"^^^ > (4-74) 
where ft^X) and 6£,„-i(A) are 

hAX) = UX + M) X I J-cl J ^^^^^ ^^^^^^^^ , (4.75) 
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The forward shift operator J-'^(A) and the backward shift operator Bi{X) are defined by 
^,(A) S M^; A + ^)-' o A.(X) o ^,(x; A) X { ; ^^CJ^^lJ:, 'd<3M . (4.77) 

def 1 H- . X f 1 : oQM, idQM 

,B,(A) = M^-, ><) ° oij,{x-X + d)x{ 1 . ,dQM , (4.78) 

and their action on the polynomials is (P£^„(x; A) =^ Pe^n{v{^)'i ^) for oQM and idQM) 

J^e{\)Pi,n{x; A) = /,,„(A)P,,„_i(x; \ + d), (4.79) 
B,(A)P,,„_i(x; A + 6) = 6,,„_i(A)P,,„(a;; A). (4.80) 

The exphcit forms of ^^(A) and B^(A) are the following: 
oQM: 

^'^'^^^'^-'l^ii-'M;^^^)), (4.81) 



idQM: 

^'^^^ ^ ^(x)Uv(x);X) (^^^^^"^ + ^i); ^ + ^y^' - ^^(^(^ - ^i); ^ + ' (4-83) 

^^(^^ - e.(,(.rA + ^) ("^^^^ ^ + ''^^^^^^^^ + ^^^^ ^^^'^ 

- y*(a;; A + id)Ce{v{x - i|); A)e-i^)(^(x), (4.84) 

rdQM: 

= r ^^ (^'^(^ + 1; A + <5) - 6(a;; A + d)e') , (4.85) 
(p{x; A + £d)^e{x + 1; A) v / 

Be{\) = ^ ^ 



B{0;X + e5)^,{x;X + S) 
X (^B(x; A + eS)^e{x; A) - D(x; A + eS)^e{x + 1; A)e-^) ^(x; A + M). (4.86) 

4.4 Second Order Equations for the New Polynomials 

The new orthogonal polynomials P^,„(77;A) satisfy second order differential or difference 
equations with rational coefficients. Here we list the explicit forms. 
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The similarity transformed Hamiltonian 'H(,{X) is defined by 

and its action on the polynomials is 

H,{X)Pt,n{x-, X) = £i,niX)Pe,nix; A). 

From fHTQil -f OOl) and f H75|) -( H76|) . we have 

Si,nW = fi,nWbi,n-i{X) = /„(A + M)6„_i(A + id) = S{n; X + id). 

The explicit forms of l-i^{X) are the following: 
oQM: 

Ht{X) = -4(02(77)^ + {c,iv,X + i6)-2c2{v)d,\ogUv; A))|^ 



- 2d2{v)ds{X, i)d, log A) - i j , 

idQM: 

, .^iiv{x + i^y,X) f ^Ar](x-t^y,X + 5) 



i,{r,{x-il);X) V i,{r,{x)-X + 6) 



e,(r^(x + zi);A) V 6(r](x);A + 5) 

rdQM: 

Ux\X) + 1; A + (5) 



^ ^ ^ ^6(^ + l;A)V 6(x;A + (5) 

+ A + ^5)M^±llA) /6(^^^liA+^ _ ^-a 
6(a:;A) V Ux;X + 6) 

For the XL and XJ polynomials in oQM we write concretely [9]: 
XLl : vdlPeAv, X) + + i + ^ - ~ 2^l^f^)d,P,Av- A) 



+ (2^^^-^ + n - l)PUm A) = 0, 
XL2 : r^dlP.^r^- X) + [g + I + \ - - 2^^||^)^.P.,n(r/; A) 

+ (-2 + n + £j A) = 0, 
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(4.95) 




2{g + \){l + ri)d,Ur^-\ + 5) 

Urn A) 

+ l{l + h-g-l)+n{n + g + h + 2^)) P^^nim X) = 0. 



(4.96) 



The zeros of the deforming polynomial ^e{rij; X) = {j = 1, . . . ,i) are the extra regular 
singularities on top of those for the original Laguerre and Jacobi polynomials. In all the four 
cases, the characteristic exponents at these extra £ zeros are the same: 



The new polynomial solutions have p = at all the extra singularities. The XJ polynomials 
provide an infinite family (n = 0, 1,2, ... ,) of global solutions of the Fuchsian differential 
equations f l4.95l) - fl4.96l) with 3 + i regular singularities. To the best of our knowledge, global 
solutions of a Fuchsian differential equation with more than four singularities had been 
utterly unknown, although in this case, the locations of the singularities are very special. 

4.5 Intertwining the Original and the Deformed Systems 

Here we present the transformation intertwining the original and the deformed Hamiltonian 
systems. It provides a simple derivation of the bi-linear expression of the new orthogonal 
polynomials in the original and the deforming polynomials. It also yields another proof of 
the shape invariance of the deformed systems. This type of transformations are sometimes 
called Darboux-Crum transformations. 

The general setting of the transformation theory is common to oQM, idQM and rdQM. 
Let us define a pair of Hamiltonians '^^^^(A) in terms of the operators ^^(A) and ^^(A)^, 



p = 0, 3. 



(4.97) 



•H^+)(A) = A(A)U,(A), •^(-^(A) ^= i,(A)A(A)t, 



(4.98) 



47 



and consider their Schrodinger equations: 

?;^f)(AMg(x;A) =4tHA)CH^;^) (^ = 0,1,2,...). (4.99) 

By definition, all the eigenfunctions must be square integrable (square summable). The pair 
of Hamiltonians are intertwined: 

•Hj+)(A)i,(A)t = Ae{X)^AiWM^y = ^(A)^H^^(A), (4.100) 
M^We'H^) = ^(A)A(A)U,(A) = H^)(A)i,(A). (4.101) 

In all the cases to be discussed here, we have Ae{X)(f)[^j^{x; A) 7^ and Ae{Xy cj)[~l^ {x; A) 7^ 0. 
This situation is called the 'broken susy' case in the supersymmetric quantum mechanics. 
This means that the two systems are exactly iso-spectral and there is one-to-one correspon- 
dence between the eigenfunctions: 

Si^j{\)=si-\X), (4.102) 
0{7j(x; A) oc A(A)0(+j(x; A), (^5+j(x; A) oc i,(A)t0(-)(x; A). (4.103) 

The strategy is to find the operators Ae and AI in such a way that "H^^^ is the original 
Hamiltonian (up to an additive constant) shown in §2.11 and ^ becomes the deformed 
Hamiltonian (up to an additive constant) demonstrated in §4.3[ Then the above formula 
(14.1031) gives the expression of the new orthogonal polynomials in terms of the original ones. 

We define operators .4.^ (A) and ^J(A) as follows: 
oQM: 

AeW = ^-d.M^;X), AeiX)^ = - d^M^; X), (4.104) 



- / A N def , ^ / / \ X N , f kx"^ + (g + £ - l)\ogx : XLl 

Mx;X) = iogUvix);X) + ^^^^^_^^^^^^_^^^^^ : XL2, XJl, XJ2 • 

(Note that ^x'^ + {g + i - 1) logx = w{ix; A + (^ - 1)<5) + const, for XLl.) 
idQM: 



(4.105) 



^(A) 'tU{eiPJv;ix; A) - e-'i^JVei^; A) ), 



AiiXy = -i{sJV,{x- A) e^P - ^F/(x; A) e"^^) , (4.106) 
y,(x; A) V{x- t(A +{£- m) ^'^f^ 7 (4-107) 
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Vnx- X) = V* {x; t(A + 1)5)) ^^M^±^^, (4.108) 

rdQM: 

^(A) = y5^(^-e^i/5^(^, i,(A)t = 1/ A) - A) e'^ (4.109) 

^ . .def6(x + liA) r -D(-(x + /? + £- 1); A + (£-1)5) :XM 

^' ^ 6(a;;A) \ D(x;t(A + (£-l)5)) : XR, XgR ^ ^ 

Then, by using the two identities of the deforming polynomials ( ]4.3p -( l474l) . ( ]4.15p -f j4.16p and 
f OSjl -f lQe]) . we can show that [Ml EHl EH [11] 

h(+)(A) = keWinX + iS + S) + /,,o(A)b,,o(A)), (4.112) 

n'--\X) = k,{X){n,{X) + /,,o(A)6,,o(A)), (4.113) 

where ^^(A) is 

dcf f 1 : XMP, XW 



^,(A) 1 : oQM, k,{X) { l^^^^^tyi \ 



- / \ \ dcf f 1 : XM, XR , , ^ , , . 

^^^^^ = I {ahd-^r^ : XgR • ^^-^^^^ 

Therefore the original system with the shifted parameters (7/(A + £<5 + <5)) and the deformed 
system (^{^(A)) are exactly isospectral. 

Based on the results (I4.112|) - (l4.113p . we have 

0(+^(x; A) = 0„(x; X + 15 + 5), ^^^^(x; A) = 0,,„(x; A), (4.115) 
4n (A) = «:,(A)(^n(A + M + 5) + /,,o(A)6,,o(A)) = /i,(A)(^£,„(A) + /,,o(A)6,,o(A)). (4.116) 

The correspondence of the pair of eigenfunctions 0^^^ (x) with their own normalisation spec- 
ified in the preceding sections is related by 



^ ^KA)0fcH^;A) ^ r 1 : oQM, idQM ^^^^ 

"^'■'^^ ' ^ v/^/,,„(A) \ y^6(l;A).,(A) :rdQM ' ' 

i,(A)t0tn(a:;A) r 1 :oQM,idQM 

VKHA)0£,n(A) I ^it{l-X)si{\) 
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By removing the effects of tlie ortliogonality weiglit functions, we define tlie operators ^^(A) 
and ^^^(A) 

. f AJ\) - f 1 : oQM, idQM 

(4.119) 

Be{X) = (l)o{x; X + £6 + 6) o —=== o i^^lx; X) x < i ■ rdOM 



(4.120) 



Tlieir explicit forms are tlie following: 
oQM: 



j-,(A) = 2[dMUv; A)|^ - d,{XMv; a + 5)) x | ; xji/m ' ^^'^^^^ 

;2^x^ -2 /C2(r/) , , / ±1 : XL1/XL2 

idQM: 

^KA) = ^ (^i(x; A + M)6(^(^ + ^i); A)e^P - t;*(x; A + iS)Uvi^ - ^i); A)e-i^) , (4.123) 

^^(^) = cf r\ 7\ (^2(a;; A + (^ - l)5)ei^ - ^;2*(x; A + (^ - 1)5)6"^^) , (4.124) 
^<.(77(x); A) V9(a;) V / 

rdQM: 

•^^(^) = , , \, f^f A + £(5)6(:r; A)e^ - <(x; A + M)6(x + 1; A)) , (4.125) 
ip{x; X + £d + 0) ^ ' 

B,{X) = \ ] U{x- A + (£ - 1)6) - v^{x- A + (£ - l)d)e-' 

Fax: X) (n T- \ 4- r — \ ]X 4- n) V 



(4.126) 

The operators ^^(A) and Bi{X) act as the forward and backward shift operators connecting 
the original polynomials P„ and the exceptional polynomials Pi^n'- 

P,{X)Pn{x- X + (.5 + 8) = knWhnix; A), (4.127) 
Bi{X)Pe,n{x; X) = 6,,„(A)P„(x; X + iS + S). (4.128) 

The former relation IKTT7\ with the explicit forms of Pi{X) iKT2T\\ . fHl23|) and f Hl25|l 
provides the explicit expressions fl4.4ip . fl4.48p and fl4.52p of the exceptional orthogonal poly- 
nomials. In terms of J^(A) and i3^(A), the relations fl4.112p - (l4.113l) become 

-B,(A) J-,(A) = niX + eS + d) + /,,o(A)6,,o(A), (4.129) 
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Pt{X)B,W = -HiiX) + /,,o(A)6,,o(A). 



(4.130) 



The other simple consequences of these relations are 

£it^W = kiWknWknW, ^n(A + id) = fe,nWknW - ^0 ( A)6,,0 ( A) . (4.131) 



Orthogonality relations fj4l3l) - fl4l4ll . (14501) - fj45B and fH755l) - fH36l) can be shown by using 
these intertwining relations [38l l37t [TT] . 

It is interesting to note that the operator ^^(A) intertwines those of the original and 
deformed systems ^(A) and ^^(A): 



AeiX + 6)AiX + eS + 6)= A(A)A(A), 
i^(A)^(A + iS + Sy = AiiXyAeiX + S). 



(4.132) 
(4.133) 



In terms of the definitions of the forward shift operators J^{X) f l2.128p . ^-'^(A) f l4.77p . ^^(A) 
( Kim , and B{X) (12:T30D . Bi{X) (Km . the above relations are rewritten as: 



s^(A + 6)PeiX + d)J^iX + £6 + 6) = s^(A) J^^A) J'^(A), 
s,(A) J-,(A)i3(A + 16 + 6) = Si{X + 6)B,{X)P,{X + 6), 



(4.134) 
(4.135) 



where s^(A) is 



Si{X) = 1 : oQM, Si{X) = ^k^{X) : idQM, s^(A) = ki{X) x { c + i-1 



c + l-l 

1 £—1 

1 — cq 



XM 
XR . 
XqR 
(4.136) 

These relations can be proven by explicit calculation with the help of the two identities of 
the deforming polynomial fjT3|) - fji^ . fl4.15p - fl4.16p and fl4.25p - fl4.26p . and provide a proof 
of fl4.73p - fl4.74p and fl4.79l) - fl4.8UI) without recourse to the shape invariance of the deformed 
system, see [3H1 ETJ E] . 

The intertwining relation offers another proof of the shape invariance of the £-th new 
orthogonal polynomials through the established shape invariance of the original polynomials 
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as depicted in the following commutative diagram: 





established 
shape invariance 



shape 
invariance 




AfiX] 




original polynomial 



new polynomial 



Two ways of proving shape invariance of the new orthogonal polynomials system. 

Historically the first members of the Xi Laguerre polynomials were discussed in the 
framework of 'conditionally exactly solvable problems' [77] in 1997. About ten years later 
the first members of the Xi Laguerre and Jacobi polynomials were constructed by Gomez- 
UUate et al [57] in 2008 in the framework of the Sturm-Liouville theory. They were rederived 
as the main part of the eigenfunctions of shape invariant quantum mechanical Hamiltonians 
by Quesne and collaborators [M]- In 2009 the present authors derived the infinitely many 
Xi Laguerre and Jacobi polynomials by deforming the Hamiltonian systems of the radial 
oscillator and the Poschl- Teller potential in terms of the eigenpolynomials of degree ^ [7]. 
The main idea was very simple. The Xi Jacobi Hamiltonian of Quesne [68] in our notation 
read 



+ 



g{g+l) h{h + l 

' ■ 2 ' 

sm X 
8{g + h + l) 



COS^ X 



{2 + g + hY 
8{2g + l){2h 



l+g + h+{g —h) cos 2x 
Its groundstate wavefunction was 

(j)ifi{x; g, h) = (sin x)^^^ (cos x) 



(sinx)^+Vcosx)'' 



{l+g+h+{g -h) cos 2x)2' 

+ g + h + {g — h) cos 2x 
l + g + h + {g — h) cos 2x 

_,p['^'-'^-'-'-'^\co.2x) 



(4.137) 



Pi' 



-1- 



cos 2x) 



(4.138) 



From this it did not take long to guess the general I form: 



<Pifl{x]g,h) 



isuixY^H cos xf+^^ 



COS 2x) 



P,. 



cos 2x) 



(4.139) 
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which was e^'^'^^'^'^ in fl4.63p for XJl. After the discovery of the second lowest members 
(£ = 2) of the XL2 family [78j, the entire XJ2 and XL2 families were constructed in [8]. 
Then the construction of the exceptional Wilson and Askey- Wilson polynomials [10] and the 
exceptional Racah and g-Racah polynomials |Tl] followed by taking shape invariance as a 
guiding principle. The Fuchsian properties and many other aspects of XL and XJ polynomials 
were demonstrated in [9]. The intertwining relations were developed in [361 133 [H]- The 
general knowledge of the solution spaces of exactly solvable (discrete) quantum mechanical 
systems governed by Crum's theorem [27] and its modifications [291 |2ll [251 [26] had been very 
helpful for the discovery of various exceptional orthogonal polynomials. Slightly different 
formulation of the Darboux-Crum transformation for the XL polynomials were reported in 
[791180]. 

Naively one expects that most exactly solvable QM Hamiltonian systems would admit 
similar shape invariant deformations. Since the 'twisting' of the parameters is essential, 
the systems corresponding to the Hermite and the Gegenbauer {g = h case of the Jacobi) 
polynomials do not admit such a deformation. So far, in addition to the examples presented 
in this review, the exceptional polynomials for the continuous Hahn [37], dual (q')-Hahn and 
little g- Jacobi [11] are explicitly known. Various exceptional orthogonal polynomials can be 
obtained from XAW and XgR polynomials by taking appropriate limits. It is a big challenge 
to construct various new orthogonal polynomials, for example, those corresponding to the 
Morse and Scarf potentials and the (q')-Hahn polynomials, and to study them in detail. 

5 Summary and Comments 

Discrete quantum mechanics, as formulated and developed by the present authors during 
the last decade, is briefly reviewed following the logical structure rather than the historical 
developments. The parallelism and contrast among the ordinary quantum mechanics (oQM) 
and the discrete quantum mechanics (dQM) are emphasised. As far as possible, one uni- 
versal formula valid for oQM and dQM with the pure imaginary shifts (idQM) and the real 
shifts (rdQM) is presented first and various ramifications follow. In one aspect, this is a re- 
formulation of the theory of bi-spectral orthogonal polynomials in the language and logic of 
quantum mechanics and every quantity is given explicitly without recourse to the moment 
problem as exemplified in ^ While it offers structural understanding of various specific 
properties of the bi-spectral polynomials, a new unified theory is established to generate all 
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the known bi-spectral polynomials from the first principle, as reviewed in ^ As evidenced 
by the discovery of various kinds of infinite families of new orthogonal polynomials reviewed 
in ^ this reformulation has been extremely successful. 

Here are a few words on how our project of discrete quantum mechanics started. The 
classical and quantum integrability were very closely related in the multi-particle dynamics 
of Calogero-Moser systems [13]. Reflecting the quantised eigenvalues etc, the corresponding 
classical quantities showed Diophantine properties at various levels [81]. In particular, the 
classical equilibrium points of the Calogero-Moser systems were described by the zeros of 
the classical polynomials, the Hermite, Laguerre and Jacobi [821 [83l [H], which in turn 
constituted the eigenpolynomials of single particle oQM. The classical equilibrium point, 
namely the minimum of the classical potential, could be better rephrased by the maximum 
of the groundstate eigenfunction (j)o{x) = e'^^^\ or the maximum of the prepotential w{x). 
Stieltjes investigated the equilibrium problems of multi-particle systems having 'logarithmic 
potential' which was nothing but the corresponding prepotential. Investigation of the discrete 
counterparts of the Calogero-Moser systems, that is, the Ruijsenaars-Schneider-van-Diejen 
[16] systems revealed similar Diophantine properties [85]. Later it was discovered that the 
classical equilibrium points of the Ruijsenaars-Schneider-van-Diejen systems were described 
by the zeros of the Wilson and Askey- Wilson polynomials [861 [87]. Then it was quite natural 
to seek for a single particle quantum mechanical formulation in which the Wilson and Askey- 
Wilson polynomials were the eigenpolynomials. Thus the first discrete quantum mechanics 
was born [12j. The factorised Hamiltonian was the guiding principle. 

The discrete quantum mechanics is still in its infancy. Although the single particle 
dynamics has been understood rather well by now, the multi-particle dynamics is virtually 
untouched. A lot of interesting problems, for example, the construction of multi-particle 
version of various results in this review, are waiting to be clarified. 

Due to the length limit, many interesting and important topics could not be covered 
in this review, for example, the limits from dQM to oQM and those of the corresponding 
poljTiomials, continuous i version of XL and XJ, etc. We refer to our papers [7]-[3H] and 
[55] and others. 
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Appendix: Symbols and Definitions 

Here are several symbols and definitions related to the (g)-hypergeometric functions, 
o shifted factorial (Pochhammer symbol) (a)„ : 

(a)n = f[{a + k-l) = a{a+l)---{a + n-l) = (A.l) 

i [a) 



k=l 

o gr-shifted factorial (g- Pochhammer symbol) (a ; q)n '■ 

n 

def 



a ; g)„ 1^(1 - = (1 - a)(l - ag) ■ ■ • (1 - ag-^). (A.2) 



k=l 

o Hypergeometric series • 



h, ■■■ ,bs 



z 



{h, ■ ■ ■ ,bs)n n\ 

It — \j 

def T-rr 



def (Qi, • • • , ar)n Z ( \ Q\ 

- 1^ (h. ... h-) ;j ' ^^-"^^ 

n=0 



where (d, • • • , ar)n = YVj=Mj)n = (ai)n • • • (Or)n- 

o g- Hypergeometric series (the basic hypergeometric series) '■ 

CO / \ ^ 

where (ai, • • • , ; g)n =^ 11^=1 (% ; ?)n = (oi ; g)n • • • (or- ; ?)n- 
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